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FURTHER  STUDIES  OF  AERODYNAMIC  LOADS 
AT  SPIN  ENTRY 


by  S.  B.  Spangler  and 
M.  R.  Mendenhall 

Nielsen  Engineering  & Research,  Inc. 


1.  INTRODUCTION 

A program  is  being  conducted  by  Nielsen  Engineering  & Research,  Inc. 
(NEAR)  for  the  Office  of  Naval  Research  to  develop  techniques  for  predict- 
ing the  aerodynamic  load  distribution  on  a fighter-bomber  type  aircraft  at 
incipient  spin  entry  flight  conditions.  The  problem  considered  is  one  in 
which  the  aircraft  is  in  the  range  of  angle  of  attack  (generally  30°  to 
45°)  in  which  flow  separation  occurs  on  the  nose  and  a steady  asymmetric 
vortex  system  is  formed.  The  vortices  induce  a side  force  on  the  nose 
and  pass  over  the  wing  and  tail  to  indue*.,  a side  force  and  rolling  moment 
on  these  components.  This  vortex  system  is  an  important  source  of 
lateral  aerodynamic  effects  causing  spin  departure  for  many  aircraft, 
which  can  occur  even  at  rero  sideslip. 

The  general  approach  to  an  analytical  description  of  the  flow  and 
the  induced  forces  and  moments  is  to  characterize  the  separated  flow 
region  over  the  nose  in  terms  of  potential  vortex  filaments  whose 
strength  and  position  over  the  nose  can  .>s  calculated.  These  filaments 
pass  over  the  wing,  body,  and  tail  and  induce  further  yawing  and  rolling 
moments.  Classical,  wing-body- tail  vortex  inteference  methods  can  be  used 
to  compute  the  positions  of  the  vortex  filaments  aft  of  the  nose  and  the 
vortex- induced  loads.  For  greater  flexibility  and  accuracy,  a vortex- 
lattice  lifting-surface  program  was  developed  which  computes  the 
asymmetric  load  distribution  inducea  on  the  wing-body  by  the  nose 
vortices. 

In  the  first  two  years  of  the  program,  the  work  has  been  primarily 
analytical.  Methods  were  developed  to  predict  the  asymmetric  vortex 
shedding  from  a circular  cross  section  nose  and  to  follow  these  vortices 
aft  over  the  wing/s trake/body/tail  to  calculate  their  interference  on 
these  aircraft  components.  The  methods  involve  a vortex'lattice  program 
capable  of  treating  asymmetric  induced  flows,  slender-body  methods  for 
tracking  the  vortices  over  the  aircraft,  and  body-tail  interference 
methods  for  calculating  vortex- induced  loads  on  the  tail  and  afterbody. 

The  vortex-lattice  program  was  subsequently  modified  to  handle  sideslip. 
This  work  is  reported  in  references  1 and  2. 
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It  was  apparent  in  applying  these  methods  to  real  aircraft  that  it 
would  be  necessary  to  consider  more  realistic  nose  shapes,  since  it  was 
clear  that  tailoring  the  shape  of  the  nose  can  improve  the  aerodynamic 
characteristics  at  high  angles  of  attack.  Since  very  little  work  lias 
been  done  on  vortex  shedding  from  noncircular  noses,  a combined  experi- 
mental and  analytical  program  was  formulated  in  the  third  year.  The 
experiments  involved  tests  in  the  V/STOL  tunnel  at  Langley  Research 
Center  and  in  the  NilAR  water  tunnel.  The  analytical  work  involved 
development  of  a vortex  shedding  analysis  for  noncircular  nose  shapes, 
including  the  effects  of  boundary  layer  separation  prediction.  This 
work  is  summarized  in  this  report.  The  report  takes  the  form  of  three 
sections  describing  the  work  and  results  and  a final  section  discussing 
and  comparing  the  results. 


2.  NASA  TEST  PROGRAM 

The  Langley  Research  Center  of  NASA  has  done  a great  deal  of  work 
on  spin  dynamics  and  aerodynamics  and.  has  produced  some  excellent  results 
on  many  facets  of  this  problem.  On  the  basis  of  their  interest,  a joint 
program  between  the  Dyn£unic  Stability  Branch  and  ONR  was  established  in 
which  an  existing  NASA  model  would  be  tested  in  the  LRC  V/STOL  tunnel  and 
simultaneous  measurements  of  the  loads  and  flow  field  would  be  obtained. 
Since  LRC  had  no  capability  to  measure  flow  fields  in  their  tunnel,  funds 
for  a probe  actuator  were  provided  for  these  tests  by  ONR.  NEAR  designed 
the  test  program  to  provide  data  suitable  for  both  the  present  ONR  con- 
tract and  NASA  use.  This  section  describes  the  tests  and  results. 

2.1  Objectives 

The  principal  objective  of  the  program  is  to  obtain  simultaneous 
measurements  of  both  the  flow  field  above  the  model  and  the  loads  on  the 
model  for  a research  configuration  having  certain  essential  features  of 
a modern  fighter-bomber  aircraft,  since  no  such  measurements  exist  (at 
least  in  the  available  literature) . Specifically,  a configuration  exhib- 
iting a noncircular  nose  and  strakes  (leading  edge  extension)  was  desired. 
The  results  were  to  be  obtained  at  low  speeds  (no  compressibility  effects) 
and  at  high  Reynolds  numbers.  The  data  were  desired  for  evaluating  a 
theory  on  vortex  shedding  from  noncircular  cross  section,  pointed  noses 
at  high  angles  of  attack.  Thus  it  was  necessary  to  measure  at  least  two 
components  of  velocity  over  a large  enough  grid  to  permit  the  vorticity 
field  to  be  well  defined, 

2.2  Model 

The  model  used  is  an  existing  NASA  model  developed  for  investigating 
high  lift  features  on  a fighter  configuration.  The  principal  character- 
istics are  shown  in  Figure  1.  The  body  has  a pointed  nose  and  a cross 
section  consisting  of  two  circular  arcs  and  straight  sides.  The  wing  is 
a clipped  delta  planform  with  a slightly  swept  forward  trailing  edge,  a 
taper  ratio  of  0.228  aspect  ratio  of  2.5,  and  a circular  arc  airfoil 
section  varying  from  6 percent  thick  at  the  root  to  4 percent  at  the  tip- 
The  model  is  geometrically  similar  to  that  described  in  reference  3 but 
approximately  2.4  times  larger. 

The  NASA  model  is  designed  to  have  canards,  which  were  not  used  in 
these  teats.  The  canard  mounts  were  used  to  support  strakes,  as  shown 
by  the  dashed  lines.  The  strakes  have  sharp  leading  edges  and  uniform 
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thickness.  An  internal  6 component  strain  gage  balance  ts  mounted  within 
the  body  to  measure  loads. 

A seven-probe  rake  was  used  to  measure  the  flow  field  above  the  body 
and  wing.  This  rake  was  mounted  on  an  actuator  located  on  the  sting 
behind  the  model,  as  shown  in  Figure  2.  On  the  hemispherical  nose  of 
each  probe  are  located  a central  hole,  four  holes  90  degrees  apart  approx- 
imately 45  degrees  off  the  axis  of  the  probe,  and  a static  hole  aft  on 
the  cylindrical  part  of  the  probe,  as  shown  below. 
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The  probes  were  individually  calibrated  to  measure  the  magnitude  and 
direction  of  the  onset  flow  up  to  flow  angles  of  60  degrees  from  the 
probe  axis  (or  alternately  three  velocity  components;  one  along  and  two 
perpendicular  to  the  probe  axis) . 

The  preliminary  design  of  the  actuator  was  done  by  NEAR,  and  the 
actuator  was  procured  by  LRC  with  ONR  funds.  The  actuator  is  an  r,  U 
type  with  a remotely  actuated  drive  in  angle  and  a manual  setting  of 
radius.  Scanivalves  were  located  on  the  sting  just  aft  of  the  actuator 
to  record  the  42  pressures  measured  by  the  probes. 

2.3  Tost  Conditions 

The  V/STOL  tunnel  is  an  atmospheric  closed-return  tunnel  with  a 
maximum  speed  capability  of  about  100  in/sec  (200  kts) . The  sting  support 
system  has  an  angle  of  attack  range  of  about  28°  and  a translation 
capability  to  place  the  model  in  the  center  (vertically)  over  the  angle 
range.  A knuckle  was  used  to  obtain  an  angle  of  attack  range  of  12°  to 
40°. 

Force  and  moment  tests  were  conducted  initially  over  the  range  from 
20°  to  40°  to  find  angles  of  attack  at  which  the  forces  wore  repeatable 
and  the  nose  separation  vortex  system  appeared  stable.  These  tests  were 
conducted  at  speeds  of  19  m/sec  and  60  m/sec.  Based  on  these  results, 
angles  of  attack  of  20°  and  35°  and  sideslip  angles  of  0°  and  10°  were 
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selected  for  the  flow  field  testing.  The  20°  case  gives  a symmetric 
vortex  pattern  with  little  side  force  and  the  35°  case  gives  an  asymmetric 
vortex  system  with  considerable  side  force.  The  flow  field  measurments 
were  all  obtained  at  speeds  of  about  60  m/sec  (200  ft/sec) . 

In  addition  to  force  and  moment  testing  and  flow  field  measurements, 
some  initial  testing  was  done  at  about  15  m/sec  with  a helium  bubble  flow 
visualisation  apparatus  to  attempt  to  locate  the  extent  of  the  vortex 
flow  field  so  the  rake  apparatus  could  be  located  to  cover  all  regions 
of  interest. 

2 . 4 Measurements 

The  measurements  consisted  of  six  components  of  forces  and  moments 
on  the  body,  some  photographs  of  the  helium  bubble  flow,  and  the  flow 
field  measurements.  The  latter  were  taken  at  3 axial  stations,  shown  on 
Figure  1.  The  first,  at  Body  Station  BS  51.9,  corresponds  to  the  inter- 
section of  the  strake  leading  edge  and  the  body.  The  second,  BS  69, 
corresponds  to  the  intersection  of  the  basic  wing  leading  edge  and  the 
body.  The  third,  BS  98,  corresponds  to  the  intersection  of  the  wing 
trailing  edge  and  the  body. 

At  the  two  forward  stations,  the  rake  was  initially  set  so  that  the 
innermost  probe  was  approximately  2.5  cm  (1  inch)  above  the  body.  The 
probe  was  then  swept  through  a ±72  degree  arc  about  the  vertical  plane. 
Readings  were  taken  every  3 degre'ss  of  arc.  The  rake  was  then  moved 
radially  outward  2.5  cm  (1  inch)  and  the  process  repeated.  Since  the 
probes  are  5 cm  (2  inches)  apart  on  the  rake,  this  procedure  provided 
data  along  arcs  every  2.5  cm  from  2.5  cm  above  the  body  to  32.5  cm  above 
the  body . 

At  the  rear  station,  the  process  was  repeated  over  a 130  degree  of 
arc.  In  addition,  the  probe  was  moved  radially  outward  to  cover  a second 
band  ranging  from  35  cm  to  6 7.5  cm  above  the  body  and  rSO  degrees  in  arc. 

It  should  be  noted  that  the  probe  measurements  turned  out  to  be 
extremely  time  consuming.  Even  though  the  V/STOb  tunnel  staff  granted 
three  additional  days  over  the  10  originally  scheduled  for  the  tests, 
less  than  half  the  data  originally  anticipated  was  actually  obtained. 

The  NEAR  Staff  engineers  received  excellent  cooperation  from  NASA  in 
these  tests. 
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2.5  Results 

The  results  of  the  tests  are  shown  in  the  figures  of  Appendix  A. 
There  are  four  sets  of  results.  They  are  the  force  and  moment  results, 
the  influence  of  the  rake  on  the  forces  and  moments,  repeatability  of 
the  flow  field  measurements,  and  vector  plots  of  the  flow  field  velocities 

The  forces  and  moments  with  no  rake  present  are  shown  in  Figures  A-1 
through  A-10.  The  normal  force  and  pitching  moment.  Figures  A-1  and  A-2, 
show  relatively  little  influence  of  Reynolds  number.  There  is  a large 
increase  in  normal  force  when  the  strakes  are  added.  The  center  of  pres- 
sure shifts  forward,  as  would  be  expected,  and  changes  with  angle  of 
attack  such  that  little  change  in  pitching  moment  occurs  with  increasing 
angle.  These  figures  also  show  the  data  from  reference  3 on  a geometri- 
cally similar  but  smaller  model  and  indicate  good  agreement  between  the 
two  sets  of  data. 

The  side  force  and  yawing  moment  results  for  various  combinations 
of  strakes  off  and  on  and  vertical  tail  off  and  on  are  shown  in  Figures 
A-3  through  A-10.  Generally  there  is  a greater  influence  of  Reynolds 
number  on  these  results,  particularly  at  the  higher  angles  of  attack. 

Since  the  viscous  effects  are  more  important  in  side  force  and  yawing 
moment,  this  is  to  be  expected.  Figures  A-3  and  A-4  show  the  smaller 
scale  results  from  reference  3 and  indicate  reasonably  good  agreement  for 
side  force  at  the  lower  angles  of  attack,  but  poor  agreement  at  high 
angles  for  side  force  and  at  all  angles  for  yawing  moment.  The  repeat- 
ability tests  made  throughout  this  angle  of  attack  range  indicated  con- 
siderable unsteadiness  in  the  forces  and  moments  (and  by  implication  the 
separation  vortex  system)  in  the  angle  of  attack  range  from  about  25  to 
32  degrees. 

The  addition  of  the  strakes  to  the  body-wing  with  no  tail  (Figures 
A-5  and  A-6)  did  not  change  the  magnitudes  of  the  force  and  moment 
appreciably  but  considerably  decreased  the  influence  of  Reynolds  number 
on  the  results. 

The  addition  of  a vertical  tail  without  strakes  (Figures  A-7  and 
A-8)  had  a large  influence  on  the  side  force  and  yawing  moment.  The 
Reynolds  number  effects  were  considerably  greater  with  the  tail  on  than 
off.  This  illustrates  the  importance  of  the  nose  vortex  system  in 
influencing  the  loads  on  the  complete  configuration. 
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Finally,  tha  addition  of  the  vertical  tail  with  strakes  on  (Figures 
A'd  and  A-10)  had  its  most  significant  effect  at  low  angles  of  attack. 

For  yawing  moment,  the  addition  of  the  vertical  tail  (compare  Figures 
A-6  and  A-10)  had  the  affect  of  reversing  tha  sign  of  the  yawing  moment 
at  angles  of  attack  above  27  degrees,  with  a large  Reynolds  number  depen- 
dency on  some  of  the  results. 

The  next  set  of  results.  Figures  A-11  and  A-12,  show  probe  inter- 
ference effects.  These  are  plots  of  the  variation  of  side  force,  yawing 
moment,  and  rolling  moment  with  rake  position  angle,  with  the  rake  inner- 
most probe  located  2.5  cm  above  the  body.  Noted  on  the  figures  are  the 
values  of  the  forces  and  moments  with  the  probe  not  present. 

For  the  forward  axial  station  (Figure  A-11) , the  probe  interference 
is  small  in  a small  angle  region  around  90  degrees  (the  probe  directly 
on  top  of  the  model) . The  interference  loads  increase  as  the  probe  is 
moved  either  direction  away  from  90  degrees  in  an  antisymmetric  fashion, 
due  to  the  lift  interference  of  the  airfoil  section  supporting  the  probes 
(see  Figure  2) . At  the  rear  station  (Figure  A-12) , the  variation  of 
loads  from  the  no-probe  values  is  about  half  that  of  the  previous  case, 
except  for  C^,  and  does  not  exhibit  the  same  antisymmetric  behavior. 

It  is  not  possible  to  assess  the  importance  of  these  interference 
loads  on  the  vortical  flow  field  measurements.  The  measurements  of 
greatest  interest  are  in  a small  angle  region  around  90  degrees  wnere 
the  interference  loads  are  the  smallest,  and  this  plus  the  repeatability 
of  the  results,  discussed  next,  lend  confidence  to  the  measurements. 

The  next  set  of  results  (Figures  A-13  and  A-14)  show  repeatability 
of  the  rake  velocity  measurements  at  two  axial  stations.  At  tha  forward 
station  (Figure  A-13) , two  sets  of  the  three  velocity  components  are 
given  for  a rake  angle  such  that  the  innermost  probes  (nos.  6 and  7)  are 
very  near  the  vortex  core.  The  agreement  is  excellent  except  for  probe 
no.  7 closest  to  the  vortex,  and  is  quite  good  at  this  point.  For  the 
rear  station  (Figure  A-14; , the  innermost  probe  is  35  cm  (14  inches) 
above  the  body.  Tne  agreement  again  is  very  good.  These  results  lend 
some  confidence  to  the  assumptions  that  the  vortex  flow  over  the  model 
is  stable  and  the  interference  effects  of  the  probe  on  the  flow  are  not 
inducing  any  large  instability  of  the  vortical  flow. 


The  final  set  of  results  (Figures  A-IS  through  A-30)  are  ”elocity 
vector  plots  prepared  by  computer  from  the  probe  data.  Each  figure  is 
identified  by  model  configuration,  axial  station,  and  flow  condition. 

The  vectors  are  drawn  to  scale  according  to  the  magnitude  of  the  com- 
ponent in  the  cross  flow  plane  (normal  to  the  model  longitudinal  axis)  , 
and  the  axial  component  is  given  by  the  number  at  the  end  of  each  vector. 
The  scale  is  dimensional,  with  the  top  of  the  body  at  4.6  inches- 
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Figures  A-15  through  A-18  show  results  at  the  most  forward  axial 
station.  For  20  degrees  angle  of  attack  (Figures  A-15  and  A-16) , the 
vortex  cores  are  very  close  to  the  body  and  are  not  well  defined.  At 
j 35  degrees  angle  of  attack  (Figures  A-17  and  A-18),  the  vortex  cores  are 
further  out  than  the  first  row  of  probe  data  and  are  reasonably  well 
defined.  One  can  note  the  asymmetry  at  zero  sideslip  (Figure  A-18)  from 
the  directions  of  the  velocity  vectors  along  the  vertical  axis.  Although 
with  the  reduction  in  the  figures,  it  is  difficult  to  read  the  axial 
velocity  magnitudes,  it  is  characteristic  of  the  vortices  at  35  degrees 
that  a high  axial  velocity  exists  in  the  core  region.  Values  as  high  as 
82  m/sec  (268  ft/sec)  occur  near  the  core,  whereas  the  typical  values 
away  from  the  core  are  49  (160)  to  55  m/sec  (180  ft/sec) . The  vortices 
tend  to  resemble  a classical  concentrated  vortex  pair. 

At  the  next  aft  station  (Figures  A-19  through  A-22) , a similar 
pattern  exists  at  20  degrees  angle  of  attack  as  occurred  at  the  more 
forward  station  at  35  degrees,  except  the  vortex  pattern  is  symnetric 
(see  Figure  A-21  in  particular) . At  35  degrees  angle  of  attack  (Figures 
A-19  and  A-20) , the  vorticity  is  more  diffuse  in  the  region  above  the 
body  and  is  quite  asymmetric. 

At  the  station  behind  the  wing,  there  are  four  sets  of  data  with 
strakes  on  (Figures  A-23  through  A-26)  and  four  without  strakes  (Figures 
A-27  through  A-30) . Unfortunately  data  were  not  obtained  for  a portion 
I of  the  outer  rake  position,  because  of  tunnel  time  limitation  and  because 
the  strake  vortex  was  more  diffuse  than  was  anticipated. 

In  general,  there  are  significant  differences  between  strake  on  and 
strake  off.  For  example,  comparison  of  Figures  A-23  and  A-27  indicates 
that  the  strake  has  a large  organizing  effect  on  the  flow.  A strong 
downwash  occurs  over  the  body  due  to  the  strake  vortices  and  there  is 
little  indication  of  the  presence  of  the  nose  vortex  system,  which  is 
probably  merged  with  the  strake  vortices.  There  is  a low  velocity  region 


20 


in  the  core  of  the  strake  vortex,  whereas  without  the  strake,  there  is  a 
large  unorganized  "dead  water"  region  over  the  wing.  The  probe  readings, 
in  fact,  in  these  low  axial  velocity  regions  are  unreliable  because  the 
angles  of  the  flow  approaching  the  probe  exceed  the  range  for  which  the 
probe  is  calibrated.  Without  the  strake,  there  is  some  indication  of  the 
presence  of  the  nose  vortex  system  at  10  degrees  sideslip.  For  example, 
in  Figure  A-28  there  is  an  indication  of  concentrated  vorticity  near  the 
vertical  axis  about  28  inches  above  the  body  axis  and  in  Figure  A-30 
there  is  a similar  indication  just  to  the  left  of  the  'vertical  axis  about 
9 inches  above  the  body  axis.  In  order  to  define  the  vorticity  field  more 
precisely,  one  would  have  to  make  detailed  vorticity  contour  integrations, 
which  have  not  been  done  at  this  point. 

Further  comments  are  made  regarding  these  data  in  Section  5. 
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3.  WATER  TUNNEL  STUDIES 

An  experimental  investigation  was  made  in  the  NEAR  water  tunnel  to 
obtain  information  on  vortex  shedding  from  noncircular  nose  shapes.  The 
program  and  results  are  described  in  this  section. 

3.1  Objectives 

The  principal  objective  of  these  tests  is  to  obtain  simultaneous 
force  and  flow  field  data  on  a series  of  noncircular  noses  characteristic 
of  real  airplane  noses.  The  nose  shapes  were  selected  to  complement 
other  parts  of  the  overall  investigation  as  well  as  to  be  representative 
of  aircraft  nose  sections.  Because  of  limited  data  acquisition  capa- 
bility, flow  visualisation  was  used  to  determine  the  vortex  locations, 
and  no  velocity  or  vortex  strength  measurements  were  made.  The  purpose 
of  the  measurements  is  to  check  other  data,  to  check  the  theory  described 
in  the  next  section,  and  to  gain  some  understanding  of  the  way  in  which 
nose  shape  affects  the  formation  and  location  of  the  vortices  shed  from 
the  nose. 

3.2  Models 

Three  models  were  fabricated  and  tested.  All  were  made  to  the  same 
base  area  (equivalent  base  diameter  of  6.25  cm  (2.5  inches)).  Two  have 
length-to-diameter  ratios  of  5 and  the  third  is  slightly  less.  The  model 
characteristics  are  shown  in  Figure  3. 

The  first  model  is  the  nose  of  the  LRC  V/STOL  Fighter  Aircraft  model 
tested  in  the  V/STOL  tunnel,  as  described  in  the  previous  section.  This 
model  was  selected  to  provide  some  data  for  comparison  with  the  LRC  data. 

The  second  model  is  a tangent  ogive  with  an  elliptical  cross  section 
having  a ma jor-to-minor  axis  ratio  of  1.44.  The  elliptical  cross  section 
is  typical  of  many  fighter  aircraft  and  is  used  on  a portion  of  the  F-5 
nose.  This  model  is  also  identical  to  one  of  the  models  tested  at  Ames 
Research  Center,  on  which  considerable  force  and  moment  data  exist 
(reference  4 ) . 

The  third  model  is  an  F-5  aircraft  nose.  '^he  length  to  equivalent 
base  diameter  is  4.4.  The  cross  section  ax-ea  shape  distribution  was 
obtained  from  Langley  Research  Center,  NASA,  and  is  identical  to  that  of 
the  model  of  reference  5.  As  a part  of  the  investigation  of  reference  5, 
a nose  alone  model  was  tested,  but  the  NASA  model  is  SMnewhat  shorter 
(in  terms  of  percent  length  of  the  overall  fuselage)  than  the  NEAR  model. 
This  modal  is  of  interest  for  two  reasons.  First,  the  aircraft  is  highly 
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spin  resistant  due  to  the  shape  of  the  nose  and  is  thus  a very  interesting 
configuration.  Secondly,  there  are  considerable  force  and  moment  data  on 
this  configuration  from  reference  5. 

All  models  were  constructed  with  a central  steel  support  tube 
fashioned  to  hold  a NEAR  3/4  inch,  5 component  strain  gage  balance. 
Aluminum  template  sections  were  located  at  a number  of  stations  along  the 
axis,  and  the  intermediate  regions  filled  with  an  epoxy  material,  with 
the  outer  surface  hand  finished  and  painted  black. 

3.3  Test  Conditions 

The  models  were  tested  in  the  water  tunnel  at  speeds  from  about  1.6 
to  6.1  m/sec  (5  to  20  ft/sec) . The  angle  of  attack  range  is  15  to  40 
degrees,  and  sideslip  angles  of  0,  i5,  and  ±10  degrees  were  tested.  Flow 
visualization  was  all  done  at  a water  speed  of  about  1.8  m/sec  (6  ft/sec). 

The  models  were  tested  in  the  0.35  by  0.5  meter  (14  by  20  inch)  test 
section  of  the  tunnel.  A photograph  of  the  arrangement  is  shown  in 
Figure  4 with  the  water  partially  drained  from  the  tunnel.  The  maximum 
blockage  was  about  7.7  percent.  No  corrections  were  made  for  blockage 
effects. 

3.4  Measurements 

Two  kinds  of  measurements  were  made;  force  and  moment  and  flow  visu- 
alization. Initially,  force  and  moment  measurements  were  made  in  the 
absence  of  any  flow  visualization  probes  or  other  disturbances.  The 
measurements  were  made  with  an  internally  mounted  three-quarter  inch 
diameter,  5 component  strain  gage  balance  (all  components  except  drag) . 

The  instrumentation  is  shown  in  Figure  4 in  the  rack  on  the  rigtit.  The 
system  is  capable  of  sampling  the  balance  output  at  about  3 samples  per 
second  (limited  by  the  printer)  and,  depending  on  the  degree  of  unsteadi- 
ness, up  to  25  samples  were  taken  and  averaged  to  obtain  a balance 
reading . 

After  the  force  and  moment  tests  were  completed,  flow  visualization 
tests  were  conducted  to  determine  the  location  of  the  vortex  cores  on  the 
leeward  side  of  the  body.  Initially  a hydrogen  bubble  technique  was 
tried.  This  involved  cementing  two  5 mil  platinum  wires  along  the  pres- 
sure side  of  the  body  about  45  degrees  around  the  curve  of  the  body  from 
the  windward  meridian.  The  wire  was  faired  into  the  body  ’,'ith  epoxy  and 
scraped  clean  over  the  location  where  bubbles  were  wanted.  This  method 
worked  only  to  show  the  separation  line  on  the  body,  because  the  bubbles 
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Figure  4.  Arrangement  of  Model  and  Instrumentation 

in  the  Tunnel. 
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would  move  along  the  body  to  the  separation  line,  accumulate  there,  and 
ultimately  break  off  to  follow  the  vortex  core.  The  bubble  density  was 
insufficient  to  show  the  vortex  cores,  and  the  method  was  limited  to 
water  speeds  of  the  order  of  0.6  m/sec  (2  ft/sec). 

The  method  used  is  one  in  which  a probe  was  placed  upstream  of  the 
model  and  used  to  supply  air  bubbles  to  the  vicinity  of  the  model.  The 
probe  consists  of  an  b'-shaped  1/8  inch  diameter  tube  entering  the  tunnel 
on  the  sidewall  upstream  of  the  test  section.  An  0.028  inch  tube  was 
soldered  to  the  large  tube  just  downstream  of  the  bend  of  the  large  tube. 
The  end  of  the  small  tube  was  soldered  and  drilled  to  0.016  inch.  By 
sliding  the  probe  through  the  fitting  at  the  tunnel  wall  and  rotating  the 
probe,  the  tip  of  the  small  tube  could  be  moved  relative  to  the  model  to 
cause  the  bubble  stream  to  impinge  the  model  just  aft  of  the  nose  tip. 

The  model  sting  support  system  was  used  to  position  the  model  approximate- 
ly correctly  relative  to  the  bubble  stream.  This  system  was  operated  at 
tunnel  speeds  up  to  about  3 m/sec  (10  ft/sec)  and  was  limited  primarily 
by  the  structural  integrity  of  the  probe. 

The  probe  was  placed  so  the  bubble  stream  impinged  on  the  model  just 
downstream  of  the  nose  tip.  The  bubbles  followed  and  showed  extremely 
well  the  vortex  core  locations.  Photographs  were  obtained  of  the  bubbles 
from  the  side  and  bottom  (leeward  model  side)  of  the  tunnel  through  the 
Plexiglas  test  section.  Examples  of  thv  photographs  are  shown  in 
Figures  5-10.  Measurements  were  made  from  the  photographs  to  obtain 
vortex  positions  at  10  positions  along  the  model  length. 

3.5  Results 

The  results  of  the  water  tunnel  tests  are  presented  in  Appendix  B. 
Included  in  these  data  are  the  basic  forces  and  moments,  an  assessment  of 
air  bubble  probe  interference,  and  vortex  position  data  on  the  four  models 
(the  elliptical  cross  section  ogive  was  tested  with  major  axis  both  hori- 
zontal and  vertical). 

The  force  and  moment  data  with  no  probe  present  are  shown  in  Figures 
B-1  through  B-16.  The  Reynolds  number  is  bt.Bed  on  the  tunnel  test  section 
velocity  and  the  equivalent  base  diameter  (the  diameter  of  a circle  having 
the  base  area) . In  some  cases  where  the  forces  were  unsteady,  a range  of 
force  meat'urad  is  shov/n.  This  occurred  only  with  side  force  and  yawing 
moment  at  the  higher  angles  of  attack.  No  corrections  were  made  for 
water  tunnel  wall  ef foots. 
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Figure  7,  Air  Bubble  Photographs  of  F-5  Nose 
at  a»40®,  S*10®. 
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Figure  9.  Air  Bubble  Photographs  of  Elliptical  Ogive 
with  Major  Axis  Horizontal  at  x = 40°,  3 * 0°. 
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The  results  for  the  LRC  V/STOL  fighter  model  nose  are  shown  in 
Figures  6-1  through  B-4.  There  is  a considerable  Reynolds  number  influ- 
ence, particularly  on  lateral  effects.  At  the  high  Reynolds  number, 
there  is  very  little  aide  force,  whereas  the  side  force  and  yawing  moment 
are  appreciable  above  30°  angle  of  attack  at  the  lower  Reynolds  number. 

The  sideslip  results  show  a linear  behavior  at  20°  angle  of  attack  and  at 
the  higher  Reynolds  number  at  35°  angle  of  attack.  At  the  lower  Reynolds 
number.  Figure  B-4  shows  a reversal  in  sign  of  the  sideslip  derivative  of 
yawing  moment.  If  the  moment  center  were  well  aft  of  the  nose  base  (as 
would  be  the  case  in  a complete  aircraft  configuration) , rather  than 
forward  of  the  base  (Figure  3),  the  reversal  would  not  occur,  but  the 
results  still  show  a significant  effect  of  the  nose  vortex  system  on  the 
center  of  pressure  of  side  force  on  the  nose. 

The  results  for  the  elliptical  ogive  with  major  axis  horizontal  are 
shown  in  Figures  B-5  through  B-8.  There  is  some  Reynolds  number  effect, 
but  not  as  pronounced  as  in  the  V/STOL  nose  case.  The  flow  becomes  very 
unsteady  at  40  degrees  angle  of  attack,  as  shown  by  the  range  of  and 
data  in  Figures  B-5  and  b-6.  The  side-force  and  yawing-moment  varia- 
tions are  reasonably  linear  with  sideslip  angle,  as  shown  in  Figures  B-7 
and  B-8. 

NASA  side-force  and  yawing-moment  data  from  reference  4 are  shown  in 
Figures  B-5  and  B-6.  These  data  show  an  opposite  ...gn  to  the  NEAR  data, 
which  is  not  unexpected,  in  that  at  zero  sideslip,  the  "hand"  of  the 
asymmetric  flow  is  determined  by  small  asymmetries  in  the  model  or  tunnel 
flow.  The  magnitudes  in  the  side-force  data  agree  reasonably  well.  Both 
sets  of  data  show  small  values  below  25°  (note  the  NEAR  data  is  plotted 
at  ten  times  the  scale  values) , but  the  NASA  data  show  higher  values  at 
the  higher  angles.  The  Reynolds  numbers  for  the  tests  are  different,  but 
there  are  insufficient  data  to  establish  the  trend  with  Reynolds  numbers. 

With  the  major  axis  vertical  (Figures  B-9  through  B-12) , there  are 
considerable  Reynolds  number  effects  at  sero  sideslip.  The  lateral 
effects  appear  at  a lower  angle  of  attack  than  was  the  case  for  the  major 
axis  horizontal,  and  the  flow  is  steadier  at  high  angles  of  attack.  The 
NASA  data  show  the  same  behavior  as  the  NEAR  data  and  are  reasonably  close 
in  magnitude,  oorisidoring  the  size  of  the  Reynolds  number  effects.  The 
variations  with  angle  of  sideslip  are  essentially  linear. 
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The  data  for  the  F-5  nose  are  shown  in  Figures  B-13  through  E-16. 

The  zero  sideslip  data  show  little  influence  of  Reynolds  number.  It  is 
remarkable  that  this  nose  configuration  shows  little  side  force  or  yawing 
moment  over  the  entire  range  of  angle  of  attack. 

Figure  B-13  shows  NASA  normal  and  side-force  data  from  reference  5. 
The  NASA  model  la  somewhat  shorter  than  tho  NEAR  model,  in  that  it  repre- 
sents 32  percent  of  the  total  fuselage  length,  whereas  the  NEAR  model  is 
41  percent  of  the  length.  Although  there  is  that  difference,  the  forces 
shown  in  Figure  B-13  agree  well  with  the  NEAR  tunnel  data.  NASA  moment 
data  are  not  presented  because  of  an  uncertainty  in  the  center  of  moments 
used,  which  could  not  be  resolved. 

The  side-force  data  in  Figure  B-15  show  a remarkable  reversal  in 
sign  at  about  34  degrees  angle  of  attack  for  both  Reynolds  numbers  tested. 
Below  34  degrees,  the  nose  force  is  away  from  the  sidewash  velocity, 
whereas  above  34  degrees,  the  side  force  is  towards  the  sidewash.  For  a 
typical  complete  aircraft  configuration  with  a center  of  gravity  aft  of 
the  nose,  the  directional  stability  also  changes  sign.  The  data  of 
Figure  B-16  do  not  show  this  effect  because  the  center  of  moments  for  the 
water  tunnel  tests  is  well  forward  of  the  model  base.  The  NASA  side  force 
data  show  the  same  trends  as  the  NEAR  data  and  the  same  crossover  near 
« - 34°. 

The  next  group  of  data  (Figures  B-17  through  B-25)  present  the 
results  of  the  flow  visualization  studies.  Examples  of  the  photographs 
obtained  with  the  air  bubbles  are  shown  in  Figures  5 through  10.  Figure  5 
shows  air  bubble  flow  over  the  V/STOL  nose.  Figures  6 and  7 show  the  P-5 
nose,  and  Figures  8-10  show  the  ogive. 

The  influence  of  the  probe  on  the  forces  and  moments  were  investiga- 
ted by  taking  concurrent  measurements  from  the  balance  with  the  photo- 
graphy. Balance  measurements  were  taken  for  the  cases  of  probe  inserted 
with  no  air  and  probe  with  air  and  compared  to  the  measurements  made  with 
the  probe  withdrawn.  The  results  of  these  measurnments  are  shown  in 
Figures  B-17  through  B-20.  Generally  these  data  show  little  influence  of 
the  probe  presence  on  the  normal  force  and  pitching  moment  (mostly  less 
than  10  percent  change)  and  a considerably  larger  effect  on  the  side  force 
and  yawing  moment.  Also  the  wake  of  tho  probe  itself  appears  to  have  a 
much  greater  influence  on  the  model  loads  than  do  the  air  bubbles,  because 
the  probe  and  probe-plus-air  interference  are  generally  very  close. 
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The  fractional  changes  in  lateral  loads  would  be  expected  to  be 
higher  than  those  in  the  vertical  plane  for  two  reasons.  First,  the 
lateral  loads  are  smaller,  particularly  at  the  low  angles  of  attack,  and 
small  differences  result  in  larger  fractional  changes.  Secondly,  the 
lateral  loads  are  caused  by  small  asymmetries  in  the  model  or  onset  flow 
and  are  much  more  sensitive  to  changes  in  the  flow  asymmetry  than  the 
vertical  plane  loads. 

On  the  basis  of  considerable  experimentation  with  the  probe  and 
model  in  the  water  tunnel  and  moving  the  probe  vertically  with  respect  to 
the  model  to  cause  the  probe  wake  and  bubbles  to  impinge  on  the  model  at 
different  locations,  it  is  felt  that  the  photographs  obtained  represent 
the  actual  vortex  system  quite  well.  In  all  cases,  the  vortex  system 
visualized  by  the  bubbles  was  quite  stable  and  did  not  jitter  or  appear 
to  be  on  the  verge  of  shifting  "hand".  In  all  cases  except  possibly  the 
elliptic  ogive  with  major  axis  vertical  at  the  highest  angle  of  attack, 
there  was  a two  vortex  system  completely  visualized  by  the  bubbles  flow- 
ing over  the  model  near  its  nose,  as  shown  in  Figure  6,  for  example,  in 
all  cases  except  the  V/STQL  nose,  the  vortex  cores  ware  well  defined,  as 
in  Figures  6 and  7,  and  the  core  positions  could  be  read  to  within  about 
■^0.10  of  the  local  body  radius.  The  vorticity  in  the  leeward  side  flow 
over  the  V/STOL  nose  appeared  more  diffuse  (Figure  5) , and  the  cores  were 
more  difficult  to  define. 

Measured  vortex  positions  for  all  oases  in  which  data  were  obtained 
are  presented  in  Table  B-1.  Selected  data  for  each  of  the  four  nose 
shapes  arS'  presented  in  Figures  B-21  through  B-25.  The  lateral  and  ver- 
tical coordinates  aro  nondimensionalized  by  the  equivalent  base  diameter. 
Also  shown  are  outlines  of  the  model  to  place  the  vortex  positions  rela- 
tive to  the  model. 

The  results  for  the  V/STOL  nose  in  Figure  B-21  are  typical  in  that 
the  vortices  lay  close  to  and  over  the  body  surface,  even  at  10  degrees 
of  sideslip.  The  lower  case  of  Figure  B-21  corresponds  to  the  photograph 
of  Figure  5,  and  there  is  some  indication  of  a second  vortex  forming  on 
the  right  aide  near  the  80  percent  station  as  the  first  vortex  bends  away 
from  the  body,  although  it  is  difficult  to  tell  precisely.  The  diffuse- 
nesB  of  the  vorticity  is  evident  by  comparing  Figures  5 and 

The  elliptic  ogive  results  of  Figures  B-23  and  b-23  indicate  a large 
vortical  height  difference  between  the  right  and  left  vortices,  even  at 


r 30  degrees  angle  of  attack/  with  the  major  axis  vettical  showing  the 

greatest  differences.  The  vortex  that  rises  the  greatest  distance  ver- 
tically also  moves  the  greatest  distance  outboard  for  zero  sideslip. 

I i The  F-5  results  shown  in  Figure  B-24  for  zero  sideslip  indicate  small 

r vertical  position  differences  and  lateral  asymmetries  for  both  30  and  40 

degrees  angle  of  attack/  which  is  consistent  with  the  small  lateral  forces 
i,  developed  on  this  configuration.  The  results  at  10  degrees  of  sideslip 

r (Figure  B-25)  bracket  the  angle  of  attack  at  which  the  side  force  direc- 

tion changes  sign.  It  is  interesting  to  note  that  for  the  smaller  angles 
of  attack/  the  vortex  system  is  sufficiently  close  to  the  body  that  when 
sideslip  is  introduced/  the  windward  vortex  position  is  strongly  influ- 
enced by  the  canopy.  It  is  possible  that  this  effect  is  responsible  for 
^ the  reversal  in  sign  of  side  force  in  going  from  30°  to  40°. 


4.  ASYMMETRIC  VORTEX  SHEDDING  ANALYSIS 

Numerous  prediction  methods  have  been  developed  in  an  attempt  to 
understand  the  separated  flow  field  and  induced  forces  associated  with 
bodies  at  high  angles  of  attack.  Generally,  these  methods  consider  a body 
of  revolution  and  represent  the  wake  as  a series  of  discrete  vortices  in 
crossflow  planes.  Typically,  empiricism  is  used  to  specify  the  strength 
and  initial  position  of  the  shed  vortices,  and  potential  flow  methods  are 
used  to  track  the  vortices  as  they  move  downstream.  One  prediction  method 
which  eliminates  the  need  for  much  of  this  type  of  empiricism  is  that  de- 
veloped by  Marshall  and  Deffenbaugh  (ref.  6) . 

In  this  reference,  a three-dimensional  steady  flow  problem  is  reduced 
to  a two-dimensional  unsteady  separated  flow  problem  for  solution.  Thus, 
the  problem  is  reduced  to  that  of  the  flow  about  a circular  body  in  a 
crossflow  plane  in  the  presence  of  discrete  vortices.  The  radius  of  the 
body  is  changing  with  time,  and  at  each  time  step,  corresponding  to  an 
interval  of  length  on  the  body,  a new  vortex  pair  is  shod  into  the  flow 
field  from  the  separation  points.  The  location  of  each  separation  point 
is  determined  by  solving  the  laminar  boundary  layer  equations  in  two 
dimensions,  and  the  strength  of  each  vortex  is  determined  by  the  vorticity 
flux  in  the  boundary  layer  at  separation  during  the  time  interval.  Some 
initial  results  from  this  method  indicate  good  prediction  of  the  normal 
forces  on  bodies  of  revolution  at  moderate  angles  of  attack.  The  analysis 
of  reference  6 is  restricted  to  the  symmetric  problem  in  which  there  are 
no  induced  side  forces. 

In  reference  7,  Deffenbaugh  and  Koerner  extend  the  prediction  method 
to  asymmetric  flow  about  bodies  of  revolution.  This  analysis  considers 
the  separation  of  turbulent  boundary  layers  through  the  use  of  Stratford's 
criterion. 

The  purpose  of  the  work  described  heroin  is  to  develop  a method, 
buued  on  the  approach  of  references  6 and  7,  to  predict  the  asymmetric 
flow  field  and  induced  forces  end  moments  on  realistic  aircraft  nose 
shapes  at  high  angles  of  attack.  The  method  is  applicable  to  nonoiroular 
cross  section  shapes,  uses  both  laminar  and  turbulent  separation  criteria 
in  a more  simplified  boundary  layer  treatment,  adds  a three-dimensional 
thickness  solution,  and  adds  consideration  of  the  axial  flow  as  vrell  as 
crossflow  component  to  improve  the  separation  prediction.  The  method  of 
analysis  and  some  predicted  results  are  presented. 
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4.1  Analysis 

4.1.1  Geornutry 

The  configurations  considered  in  this  study  of  vortex  shedding  are 
bodies  alone;  that  is,  bodies  without  fins  or  control  surfaces.  Thus, 
the  method  is  suitable  for  use  on  long  slender  missile  configurations  as 
well  as  shorter  aircraft  noses.  The  bodies  need  not  be  axisymmetric  as 
it  is  only  necessary  that  the  coordinates  describing  the  cross-sectional 
shape  be  defined  at  a number  of  axial  stations.  Neither  top-bottom  or 
right-left  symmetry  is  required,  but  there  are  some  shape  restrictions 
Imposed  by  the  numerical  methods  used.  Since  a Fourier  series  representa- 
tion of  the  transformation  function  is  used,  sharp  corners  and  long  flat 
body  panels  in  the  cross  section  are  to  bo  avoided.  Discontinuous  axial 
distributions  of  cross-sectional  areas  should  also  be  avoided  because  of 
resulting  numerical  difficulties. 


4.1.2  Noncircular  bodies 


The  solution  for  the  two  dimensional  potential  flow  around  a circle 
in  the  presence  of  a uniform  flow  and  external  vortices  is  well  known  and 
is  documented  numerous  places  in  the  literature,  for  example,  references 
6 and  8.  VJhen  the  body  cross-sectional  shape  is  noncircular,  the  poten- 
tial solution  is  not  readily  obtainable.  The  procedure  used  to  handle 
noncircular  shapes  is  to  find  a conformal  transformation  which  will  map 
every  point  on  or  outside  the  noncircular  body  to  a corresponding  point 
on  or  outside  a circular  body.  Thus,  the  potential  flow  solution  can  be 
written  for  the  circular  body  and  transformed  to  the  noncircular  body. 


For  some  simple  shapes  like  an  ellipse,  the  transformation  to  the 
circle  is  relatively  simple  and  can  be  carried  out  analytically.  As  the 
noncircular  shape  becomes  more  complex,  the  transformation  cannot  be 
done  analytically  and  some  numerical  procedure  must  be  used.  The  numer- 
ical procedure  chosen  is  a form  of  the  Theodorsen's  transformation  given 
by  the  relation 
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B8  illustrated  in  figure  11.  For  each  point  a on  the  real  body  cross 
section  there  is  a corresponding  point  v on  the  circle. 

The  Fourier  coefficients,  and  in  equation  (1)  are  obtftined 
from  a numerical  Fourier  analysis  procedure  developed  by  A.  Jameson  for 
use  in  transforming  arbitrary  airfoil  shapes  to  a circle  (ref.  9).  Given 
the  Fourier  coefficients,  the  mapping  of  known  points  in  one  plane  to 
corresponding  points  In  the  second  plane  is  carried  out  numerically  as 
described  in  Appendix  C. 

4.1.3  Equations  of  motion 

In  this  section,  the  equations  of  motion  of  a shed  nose  vortex  in  ^he 
presence  of  an  arbitrary  number  of  other  free  vortices  in  the  vicinity  of 
a noncircular  body  are  derived.  The  transformation  of  the  nonoircular 
body  to  the  circle  xu  given  by  equation  (1) , and  the  mapping  of  corres- 
ponding points  from  one  plane  to  the  other  is  described  in  Appendix  C. 

In  the  circle  (v)  plane,  let  the  position  of  a vortex,  r^,  be  denoted 


n n 


and  the  image  of  V is  located  at 

n 


The  complex  potential  in  the  circle  plane  is 


W(v)  “ 0 + iij< 


and  the  velocity  field  in  the  circle  plane  is 


v ' - iw ' « 


dW(v) 


(6) 


Real  Plane 


Figure  11.  Body  cross  Section  Nomenclature 
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The  complex  potential  of  the  flow  in  the  circle  plane  is  described  by 
the  following  components s 

(1)  crossflow  due  to  a, 

(v)  ■ - iv  V sin  Cl  (10) 

(2)  crossflow  due  to  3, 

Wg(v)  = - V V sin  B (11) 

(3)  cylinder  in  a-flow  (doublet), 


W (v)  = i —1  V sin  a 

(4)  cylinder  in  S-flow  (doublet) , 


W (V)  - - V sin  B 

B 


(5)  expanding  cylinder, 


W*  (v)  - r,  (In  v)  V COB  a 


(6)  vortex  outside  cylinder. 


(7)  image  of 
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The  total  potential  for  the  flow  in  the  circle  plane  is  the  sum  of 
equations  (10)  to  (16)  and  may  be  written  as 
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where  N ■ total  number  of  shed  vortices. 

The  complex  potential  governing  the  motion  of  vortex  in 

the  real  (o)  plane  can  be  written  as 
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v;here  the  last  term  represents  the  potential  of  in  the  o-plane 

which  has  been  removed  from  the  potential  to  eliminate  the  singularity  at 

m 


The  velocity  of  is 
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Note  that  the  last  two  terms  in  equation  (18)  are  singular  at  ° 

V “ Combining  the  last  two  terms  of  (18)  and  talcing  the  limit  results 

in  the  following. 
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The  details  of  this  derivation  are  carried  out  in  Appendix  A of  reference 

8. 


From  (20), 
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and,  from  (1), 
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where  L is  the  number  of  Fourier  coefficients  required  to  describe  the 
transformation.  Therefore,  from  (21)  and  (22) 
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The  velocity  at  vortex  is  found  in  the  following  manner. 
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We  will  assume  the  following  form  for  equation  (24) i 
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Each  term  in  (25)  represents  a velocity  component  in  the  a plane.  The 
first  term  represents  the  flow  about  a body  in  an  a-type  crossflow, 
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The  second  term  represents  the  flow  about  a body  in  a 8-type  crossflow 
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The  third  term  is  due  to  the  change  in  body  radius  and  is 
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The  next  term  represents  the  influence  of  all  vortices  and  their  images, 
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The  next  term  is  due  to  the  image  of  r . 
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The  last  term  in  (23)  represents  the  last  two  terms  of  (18)  and  Is  writ- 
ten, with  (20)  and  (23)  as 
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Thus,  equations  (25)  through  (31)  specify  the  complete  velocity  field  at 
vortex  r„  in  the  real  plane. 
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The  differential  equations  of  motion  for 
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Therefore,  the  two  equations  which  must  be  integrated  along  the  body 
length  to  determine  the  trajectory  of  are 
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There  are  a pair  of  equations  like  (34)  and  (35)  for  each  vortex  in  the 
field.  Thus,  as  new  vortices  are  shed,  the  number  of  equations  to  solve 
increases . 

The  above  differential  equations  of  motion  for  the  vortices  are 
solved  numerically  using  a method  which  automatically  adjusts  the  step 
size  to  provide  the  specified  accuracy. 

4.1.4  Surface  pressure  distribution 

The  surface  pressure'  distribution  on  the  body  is  required  for  two 
purposes!  to  calculate  the  distribution  of  forces  on  the  body  and  to  cal- 
culate the  separation  points  on  the  body.  The  surface  pressure  coefficient 
is  determined  from  the  unsteady  Bernoulli  equation  as  follows. 
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Defining  the  pressure  coefficient  as 
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(37) 
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aa  the  total  local  velocity  (including  V) > equation  (36)  becotnea 
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A difficulty  ariaaa  in  handling  the  unsteady  term  on  the  left  side  of 
(36).  In  reference  6,  the  pressure  coefficient  was  determined  within  an 
un)(nown  constant  as  a result  of  this  term^  and  a )«nown  pressure  coeffici- 
ent at  one  point  on  the  body  was  used  to  evaluate  the  constant.  In  the 
present  method*  the  two-dimensional  doublet  in  the  crossflow  plana  repre- 
senting the  expanding  radius  is  replaced  with  a three-dimensional  source- 
sink  distribution  representing  the  entire  body  (ref.  10) . This  modified 
velocity  potential  approaches  zero  at  Infinity*  thus  allowing  calculation 
of  an  absolute  pressure  coefficient.  The  details  of  this  derivation  are 
presented  in  Appendix  D. 

The  unsteady  term  in  equation  (36)  is  written  as  follows. 
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where  is  the  dimensionless  source  strength  and  J is  the  total 

ber  of  sources  describing  the  body. 


(39) 
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The  unsteady  term  of  equation  (39)  is  written  in  terms  of  the  trans- 
formed body;  that  Is,  the  axisymmetrlc  body.  This  is  necessary  because  of 
the  use  of  a simple  source-sink  distribution  on  the  axis  to  model  the  body 
growth.  We  will  assume  that  the  unsteady  term  is  unchanged  between  cor- 
responding points  on  the  circle  and  the  real  body.  The  velocity  compon- 
ents used  in  equation  (38)  are  the  actur.l  surface  velocities  on  the  real 
body. 

4.1.5  Boundary-layer  separation 

The  approach  taken  in  this  analysis  requires  that  the  separation 
points  on  the  body  at  each  cross  section  be  known  in  order  to  specify  the 
strength  and  initial  position  of  the  uhed  vortices.  In  reference  6,  the 
laminar  boundary- layer  equations  are  solved  to  determine  the  separation 
points.  In  the  present  analysis,  separation  criteria  based  on  the  poten- 
tial pressure  distribution  are  utilized  to  predict  both  laminar  and  tur- 
bulent separation. 

uaminar  separation— A separation  criterion  proposed  by  Stratford  in 
reference  11  states  that  the  laminar  boundary  layer  separates  when 

’ 0.087  (40) 

This  criterion,  with  a slightly  different  value  of  the  constant  (0.102), 
is  discussed  in  reference  12.  For  purposes  of  this  analysis,  the  original 
constant  shown  in  equation  (40)  will  be  used.  Also,  C is  the  run  length 
of  the  laminar  boundary  layer  measured  from  a false  origin. 

To  apply  equation  (40)  on  the  body  surface  in  the  crossflow  plane, 
some  redefinition  of  the  pressure  coefficient  is  required.  Stratford's 
criterion  is  based  on  a flow  in  which  the  boundary  layer  develops  on  a 
flat  plate  in  a constant  pressure  region  followed  by  an  adverse  pressure 
region  in  which  the  boundary  layer  separates.  This  is  illustrated  in  the 
following  sketch. 


The  point  at  which  the  local  pressure  starts  to  increase  is  denoted  as 
which  corresponds  to  the  minimum  pressure  point  on  a body  in  a crossflow 
plane.  Defining  a new  surface  pressure  coefficient,  C which  must  be 
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zero  at  and  Inoreaaea  downstream,  the  surface  pressure  coefficient 

required  In  equal;lon  (40)  can  be  written  as 

Cp'  . 1 - Ml) 

where  Is  the  potential  velocity  at  the  minimum  pressure  point. 

The  specification  of  the  location  of  the  false  origin  for  a laminar 
boundary  layer  on  a circular  body  in  a crossflow  is  determined  as  follows. 
As  shown  in  the  sketch  below 


/ "*0 

V 


♦v  sin  a. 


the  minimum  pressure  point  is  denoted  m and  the  false  origin  is  a dis- 
tance Zq  upstream  of  m.  From  figure  12.11  of  reference  13,  the  dis- 
placement thickness  of  the  laminar  boundary  layer  at  m,  given  by  the  von 
Karman-Pohlhaustui  solution,  is 


where  I is  the  diameter  of  the  circle.  The  displacement  thickness  on  a 
flat  plate  boundary  layer  with  constant  pressure  equal  to  the  minimum 
pressure  at  m on  the  cylinder  is 

1.721  Zn 


where  Zq  the  flat  plate  run  length.  Assuming  that  6*  ftom  equa- 
tions (42)  and  (43)  are  equal  and  that  is  the  distar 

of  the  cylinder  between  m and  the  false  origin,  we  get 

Co  - 0.49  rg 


is  the  distance  on  the  surface 
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ThuBf  tha  run  length/  C/  in  equation  (40)  is  measured  from  a point  which 
is  a length  Kq  upstream  of  the  minimum  pressure  point  on  the  body. 

It  can  be  shown  by  a similar  analysis  that  the  false  origin  on  an 
elliptic  cross-sectional  body  is  located  about  the  same  length  upstream 
of  the  minimum  pressure  point  as  on  a circular  body.  It  will  also  be 
shown  later/  that  the  location  of  the  false  origin  is  not  a crlbioal  fac- 
tor in  the  prediction  of  laminar  bound Ary- layer  separation  using  equation 
(40) > therefore/  the  length  given  by  equation  (44)  will  bo  used  to  specify 
the  false  origin  on  all  shapes  considered  in  the  current  work. 


Turbulent  separation — separation  oritorion  for  turbulent  boundary 
layers  is  proposed  by  Stratford  in  reference  11  which  states  that  separa- 
tion occurs  when 
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The  false  origin  is  specified  by 


(45) 
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assuming  a fully  turbulent  boundary  layer.  In  equation  (47)  , u^  is  the 
local  surface  velocity  and  the  integration  is  carried  out  from  the  stagna- 
tion point  to  the  minimum  pressure  point.  In  equation  (45) , is  the 
modified  pressure  coefficient/  C^/  from  equation  (41)  , and  as  in  the  lam- 
inar criterion/  ^ is  measured  from  the  false  origin. 


Three  Dimsnsionality  Rffeots — The  Stratford  criteria  are  based  on 
two  dimensional  data  and  wor}t  well  for  circular  bodies  in  orossflowr  i.e. 
90°  angle  of  attack.  For  inollned  bodies  at  angles  of  attack  leas  than 
90°/  an  axial  flow  component  i*  added  to  the  crossflow  component/  the 
boundary  layer  run  lengths  are  longer/  and  separation  tends  to  occur  a 
shorter  distance  beyond  the  minimum  pressure  points  nn  the  body.  An 
approximate  oorreotion  for  this  effect  is  obtained  by  multiplying  the 
right  hand  aides  of  equations  (40)  and  (45)  by  sin  i.e./  for  the 

turbulent  case 


P(0  " 0.35  sin 


(48) 
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4.1.6  Separated  wake 

The  potential  model  of  the  separated  wake  is  made  up  of  a large  num- 
ber of  individual  vortices,  one  shed  from  each  separation  point  on  either 
aide  of  the  body  at  each  time  step.  The  details  of  the  specification  of 
each  individual  vortex  is  described  in  the  following  sections. 


Vortioitv  flux— In  the  same  manner  as  in  reference  6, 
flux  across  the  boundary  layer  at  the  separation  point  is 


the  vorticity 


(49) 


assuming  no  slip  at  the  wall.  In  a given  time  step  ^t,  which  is  related 
to  the  distance  between  shedding  points  by  the  relationship 
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the  vorticity  flux  is  summed  into  a single  point  vortex  whose  strength  is 
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At 


(51) 


Since  the  surface  velocity,  U,  is  dependent  on  position  on  the  body, 
the  strength  of  the  shed  vortex  depends  on  the  location  of  the  separation 
point  on  each  side  of  the  body.  The  strength  of  a vortex,  once  determined 
from  equation  (51) , is  not  changed  at  any  point  downstream;  however,  the 
net  vorticity  in  the  flow  field  is  constantly  changing  because  new  vor- 
tices are  shed  at  each  Interval  of  time.  For  the  present  work,  the  inter- 
vale between  shed  vortices  are  constant  over  the  entire  length  of  the 
body.  The  oalculution  procedure  sheds  a vortex  from  each  side  of  the  body 
at  each  time  step  if  a separation  point  is  found  by  the  Stratford  criter- 
ion  described  in  a previous  section. 

In  the  calculation  procedure,  the  vortex  strength  is  referenced  to 
the  free-stream  velocity  and  equation  (31)  is  written  as 


r 

7 


1 ful 
I W 


(VAt) 


(52) 


Initial  position — As  the  vortices  are  shed  into  the  outer  flow,  they 
muct  be  given  some  initial  position  from  which  they  begin  their  motion. 

In  reference  6,  the  vortices  are  placed  at 
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Which  is  illustrated  in  the  sketch  below. 
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The  vortex  and  Ita  image  are  placed  such  that  the  surface  velocity  U ia 
exactly  eanoellad  by  the  induced  effect  of  the  vortices.  Equation  (53)  is 
based  on  the  assumption  that  mj^<<rQ,  which  is  why  the  vortex  and  its 
image  are  equidistant  from  the  body  surface  the  curvature 

of  the  surface  is  considered,  the  initial  vortex  position  outside  a circu- 
lar body  is  given  by 


2rQ(r/V) 

"'k  ■ Jirr^  (U/V)  - (r/V) 


and  the  image  position  is  specified  by 
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The  positions  given  by  equations  (54)  and  (55)  also  satisfy  the  condition 
that  the  net  velocity  on  the  surface  at  the  separation  point  is  xero. 

It  can  be  seen  from  either  equation  (53)  or  (54)  that  the  initial 
vortex  position  is  a function  of  the  location  of  the  separation  point  and 
the  length  of  the  time  step,  in  the  course  of  this  investigation,  it  was 
found  that  a reasonable  length  time  step  which  would  produce  good  detail 
in  the  wake  would  also  cause  the  initial  vortex  positions  to  be  very  close 
to  the  body  surface.  The  effect  was  that  the  vortex  was  "captured"  by  its 
image  and  did  not  move  out  into  the  flow  field  away  from  the  body.  This 
difficulty  was  overcome  by  arbitrarily  placing  the  vortex  at  a point 
specified  by 


«"k  • fo 


(56) 


whenever  from  equation  (54)  was  less  than  5 percent  of  the  radius. 

Vortex  diffusion— The  velocity  induced  by  a point  vortex  in  a cross- 
flow plane  can  be  written  as 
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(57) 
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where  ■ 1 for  a potential  vortex.  The  major  difficulty  with  equation 
(57)  oooure  when  t,  the  dietanoe  from  the  vortex  to  the  point  of  interest* 
is  small  and  the  induced  velocity  is  large.  In  the  limit  of  trying  to  cal- 
culate the  induced  velocity  at  the  center  of  the  vortex  (r-0) * the  result 
is  singular.  The  induced  velocity  from  a potential  vortex  is  illustrated 
in  figure  12. 


Real  vortices  do  not  exhibit  this  singularity  effect  but  tend  to  dif- 
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A representation  of  this  effect  it  to  let 
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V,  ■>  effective  jtinematio  viscosity 

Equation  (58)  includes  a viscous  effect  through  the  unit  Reynolds  number 
and  an  aging  effect  through  the  ^x  term.  The  Reynolds  number  has  been 
modified  by  use  of  an  affective  kinematic  viscosity  to  permit  larger 

viscosities  to  be  used  to  increase  the  diffusion  effect.  In  all  the  cal- 
culations made  to  date*  however*  this  feature  was  not  used*  and  the  same 
unit  Reynolds  number  was  applied  to  equation  (58)  as  was  applied  to  cal- 
culationsof  the  boundary  layer  charaateristics.  The  induced  velocity  from 
equation  (57)  with  the  Pq  factor  included  is  shown  in  figure  12. 


Vortex  reduction  factor — It  was  shown  in  reference  6 that  better 
agreement  between  experiment  and  theory  would  be  obtained  if  only  a frac- 
tion of  the  vorticity  produced  by  boundary- layer  separation  is  retained  in 
the  wake.  The  fraction  that  appears  to  give  best  agreement  is  0.6.  In  the 
present  investigation*  the  vorticity  shedding  factor  is  fixed  at  0.6  unless 
otherwise  noted.  This  factor  is  used  in  the  right-hand  side  of  equation  (52) . 


Vortex  coalescence — ^As  the  calculation  proceeds*  the  number  of  vor- 
tices in  the  wake  steadily  increases*  with  an  associated  increase  in  com- 
putation time  and  costs.  In  an  effort  to  decrease  the  computation  time 
without  saorif icing  the  accuracy  of  the  predictions,  certain  pairs  of  vor- 
tices are  combined  into  one  vortex  to  reduce  the  total  number  of  vortices 
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Figure  12.  Viscous  Vortex  Model. 
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in  the  field.  The  coinblnation  of  two  vortices  is  carried  out  as  follows. 
The  strength  of  the  resulting  vortex  is 
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and  the  combined  vortex  is  located  at 


'1+2 


^1 

^1 

^2 

^2 

^1 

^2 

(61) 


■i+2 


^1 

1^2 

"2 

'1 

+ 

^2 

(62) 


The  objective  is  to  reduce  the  number  of  vortices  without  malting  an  appre- 
ciable change  In  the  flow  field;  in  particular,  without  cninging  the  in- 
duced flow  field  near  the  body  surface. 

In  keeping  with  this  objective,  certain  rules  are  oonsiderea  in  com* 
bining  vcttices.  Only  vortices  of  lllte  sign  are  combined,  and  some  effort 
is  made  to  combine  vortices  which  have  about  the  same  age.  The  vortices 
to  be  combined  should  be  in  close  proximity.  A rule  of  thumb  is  that  vor- 
tices that  are  closer  than  the  two  most  recently  shed  vortices  are  candi- 
dates for  combination.  Generally,  combination  of  vortioes  is  not  consid- 
ered until  a large  field  of  vorticity  is  developed,  and  then  a second  com- 
bination is  not  carried  out  until  some  point  far  downstream  where  another 
large  field  of  vorticity  has  developed. 

Combination  of  vortices  has  not  been  restricted  to  vortex  pairs,  and 
if  three  or  four  vortioes  are  wibhln  bhe  specified  spacing,  they  are  com- 
bined. In  the  present  form  of  the  program,  the  combination  of  adjacent 
vortices  is  done  automatically  with  logic  ir.  the  program.  The  user  of  the 
program  may  select  the  distance  within  which  the  two  (or  more)  vortices 
must  lie  to  be  combined,  and  this  distance  is  input  to  the  program. 

4.1.7  Forces  and  moments 

The  forces  and  moments  on  the  body  can  be  computed  using  two  different 
procedures.  The  first  is  an  integration  of  the  pressure  distribution 
around  the  body,  and  the  second  in  a combination  of  slendar-body  theory  and 
the  vortex-impuise  theorem.  Since  the  pressures  are  calculated  for  pu’  • 
poses  of  determining  the  location  of  separation,  the  pressure  integration 
method  13  used. 
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At  a specified  station  on  the  body,  the  normal-fcrce  coefficient  on 
a Ax  length  of  the  body  is 
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The  total  normc ) force  on  the  body  is 
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The  pitching-motiiu/jt  coefficent  is 
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and  the  center  of  pressure  Is  located  at 
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Similarly,  the  side  force  on  an  element  of  length  Ax  is 

f2it 


„ AS  -1 
“ qZ  - T 


J 


Cp  rg  sin  3 dS 


(66) 


(67) 


The  total  aide  force  is 
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The  yawing-moment  coefficient  is 
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and  the  center  of  pressure  of  the  side  force  is  located  at 
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4.2  Results 

A number  of  results  were  obtained  on  circular  or  near>circular  bodies 
to  check  the  methods.  These  include  two-dimensional  and  three-dimensional 
configurations.  Typical  results  are  described  in  thiu  section. 

4.2.1  Two-dimensional  cases 

The  discussion  in  the  previous  section  pointed  out  the  dependence  of 
the  predicted  results  on  the  separation  location.  Since  the  predicted 
separation  points  are  obtained  from  two  different  criterion  by  Stratford 
for  either  laminar  or  turbulent  flow,  the  validity  of  these  criterion  and 
the  manner  in  which  they  are  applied  were  investigated. 

The  Stratford  laminar  separation  criterion  was  applied  to  the  poten- 
tial pressure  distributions  on  two,  two-dimensional  cylinders.  The  first 
cylinder  is  circular  in  cross  section  and  the  second  is  a 2il  ellipse  in 
cross  section.  Applying  the  Startford  criterion  from  the  minimum  pressure 
point  without  including  a false  origin  results  in  separation  points  as 
noted  on  Figure  13.  When  the  false  origin  is  included,  the  predicted  sep- 
aration points  move  upstream  approximately  2*^  of  arc  on  the  circular  cyl- 
inder (Figure  13a]  and  approximately  3^  of  arc  on  the  ellipse  (Figure  13b) 
The  computed  separation  points  taken  from  reference  13  (Figure  12.11}  are 
about  2°  of  arc  upstream  of  these  latter  Stratford  results. 

Next,  Stratford's  turbulent  separation  criterion  was  applied  to  a two 
dimensional  circular  cylinder  from  reference  14.  A measured  pressure  dis- 
tribution indicated  separated  flow  on  the  back  side  of  the  cylinder.  The 
measured  pressure  distribution  was  extrapolated  through  the  separated  re- 
gion to  approximate  a potential  preusure  distribution  for  use  in  the  Strat 
ford  criterion.  The  criterion  indloated  separation  (F^  > 0.40)  at  0 > 100 
As  noted  in  reference  14,  saparation  was  predicted  to  occur  at  @ 101° 

and  separation  was  observed  at  S "=  110°. 

4.2.2  Three-Dimensional  cases 

The  Stratford  turbulent  separation  criterion  was  applied  to  an  ogive- 
cylinder  at  a u 20°  (ref.  15).  Using  the  measured  pressure  distriLotion 
at  an  axial  station  on  the  cylinder  6 diameters  aft  of  the  nose,  which 
indicated  separation  at  0 «=  120°,  the  Stratford  factor  was  computed 
through  the  separation  point.  The  value  at  separation  was  “ 0.127 
which  is  approximately  equal  to  0.37  sin  Og.  This  result  tends  to  verify 
the  suggested  modifioation  of  the  separation  factor  by  the  sine  of  the 
angle  of  attack,  as  shown  in  equation  (48) . 
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Figure  13.  Comparison  of  Stratford's  Laminar  Separation  Criterion 
with  Approximate  Solution  of  Laminar  Boundary  Layer  Equation. 

b)  2tl  Ellipse 
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The  prediction  method  was  applied  to  the  three-dimensional  case  of 
the  ogive-cylinder  of  reference  15  to  evaluate  pressure  distributions, 
vortex  positions,  and  loads.  The  configuration  has  a circular  cross  sec- 
tion and  a three-caliber  ogive  nose  followed  by  a 7.7-caliber  cylindrical 
afterbody.  The  calculations  were  made  for  an  angle  of  attack  of  20°  and 
a free  stream  Reynolds  number  based  on  diameter  of  4.4  x 10^. 

The  predicted  pressure  distributions  at  two  axial  stations  on  the 
body  are  compared  with  measured  pressure  distributions  in  Figure  14.  The 
agreement  near  the  nose  is  excellent.  It  should  be  noted  that,  in  this 
region,  the  "unsteady"  portion  of  the  predicted  pressure  coefficient  has 
a large  effect  on  the  predicted  pressure  distribution;  therefore,  the  good 
agreement  exhibited  near  the  nose  indicates  that  the  chosen  approach  for 
predicting  the  unsteady  part  is  reasonable.  The  agreement  between  measured 
and  predicted  pressure  distributions  on  the  cylindrical  section  of  the  body 
is  not  as  good  as  nearer  the  nose.  The  irregular  behavior  in  the  predicted 
curve  near  3 ■ 160°  is  caused  by  interference  from  shed  vortices.  The 
lack  of  a corresponding  dip  in  the  data  may  indicate  that  the  shed  vortices 
are  not  allowed  to  diffuse  at  the  correct  rate.  The  predicted  results  were 
obtained  assuming  laminar  separation  at  F„  - 0.025. 

The  comparison  on  vortex  position  is  shown  in  Figure  15  at  a station 
7.5  diameters  aft  of  the  nose.  The  right  side  of  the  figure  shows  the 
positions  of  the  individual  vortices  shed  from  the  separation  line  along 
the  body.  The  left  side  shows  the  predicted  center  of  vorticity,  together 
with  a measured  value.  The  two  are  in  reasonably  good  agreement. 

Figure  16  shows  the  results  for  section  normal  force  coefficient. 

There  is  little  vortex-induced  effect  on  the  ogive  nose,  and  the  load  in 
this  region  is  primarily  due  to  the  expanding  nose  cross  section.  The 
load  distribution  on  the  cylindrical  portion  of  the  body  is  due  entirely 
to  separation.  The  predicted  values  are  generally  low  except  for  the 
first  2 calibers  of  length  behind  the  shoulder. 

Finally,  the  prediction  method  was  applied  to  an  ogive-cylinder  for 
which  laser  velocimeter  data  are  available  to  compare  velocity  fields  in 
the  separated  region  over  the  body.  The  tests  were  conducted  in  the  Ames 
Researchi  Center  7-  by  10-foot  low  speed  tunnel  on  a six  inch  diameter  model 
having  a 3.5  caliber  ogive  nose  and  a 3.5  caliber  cylindrical  afterbody. 

The  freestream  Reynolds  number  based  on  diameter  was  1.9  * 10  . The  test 
conditions  and  results  are  described  in  reference  16. 
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Figure  15.  Predicted  and  Measure'll  Vortex  Positions  on  an 
Ogive-Cylinder  at  x/doo  “ 7.5. 


Figure  16.  Predicted  and  Measured  Section  Normal 
Force  Coefficient  on  an  Ogive-Cylinder. 


The  theory  was  applied  to  this  case  using  the  turbulent  separation 
criterion  and  the  sin  a factor  as  given  in  equation  (48).  In  order  to 
introduce  an  asynunetrio  disturbance  into  the  calculation,  the  separation 
points  on  the  flanks  of  the  body  over  the  initial  12  percent  of  the  model 
length  were  shifted  2 degrees.  Downstream  of  the  12  percent  length  sta- 
tion, the  separation  points  as  calculated  were  used. 

Results  for  the  downwash  distribution  at  the  70  percent  axial  station 
are  presented  in  Figures  17  and  18  . The  angle  of  attack  is  37.5  degrees, 
Figure  17  shows  predicted  and  measured  vei.ocities  along  a horizontal  line 
slightly  above  the  body  surface  and  Figure  18  shows  results  along  the 
vertical  line  passing  through  the  body  center.  The  theory  in  both  cases 
agrees  quite  well  with  the  data.  The  vortex  system  is  clearly  asymmetric 
at  this  station,  and  the  theory  shows  the  asymmetric  features. 
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Pigura  17.  Predicted  and  MeaBured  Downwash  Oiatrlbution  Along  a 
Horlaontal  Line  Juat  Above  the  Body  of  an  Oglve-Cyllndor 

at  Low  Speed. 
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5.  DISCUSSION  OF  RESULTS 

Several  areas  of  the  investigation  des«;ribed  in  previous  sections 
are  worthy  of  discussion • These  have  to  do  with  the  quantitative  and 
qualitative  observations  of  the  water  tunnel  tests > .comparisons  between 

I 

theory  and  various  seta  of  data  made  possible  by  the  coordinated  tests, 
and  observations  on  the  status  of  the  nose  vortex  shedding  analysis  made 
possible  by  the  comparison. 

5.1  Nose  Shape 

One  shape  feature  that  is  suggested  by  the  water-tunnel  data  as 
being  important  for  tailoring  lateral  loads  on  aircraft  noses  at  high 
angles  of  attaok  is  that  which  affects  the  height  of  the  vortices  above 
the  body.  The  F-5  nose  has  a vortex  pair  that  lays  relatively  close  to 
the  body  at  angles  of  attack  up  to  40  degrees.  If  this  is  the  case,  then 
as  the  nose  is  sideslipped  and  the  vortices  move  laterally  with  respect 
to  the  body,  the  windward  vortex  is  over  the  center  of  the  body  and  can 
be  affected  by  the  presence  of  the  cockpit  canopy.  The  influence  could 
take  the  form  of  lateral  and  vertical  deflection  of  the  core  and  possibly 
a change  in  axial  pressure  gradient  which  could  change  the  vortioity  dis- 
tribution in  the  crossflow  plane.  The  latter  is  suggested  by  the  differ- 
ence in  appearance  of  the  vortioes  between  Figures  6 and  7 , which  appear 
to  be  quite  concentrated  at  30  degrees  angle  of  attack  and  more  diffuse 
over  the  canopy  area  at  40  degrees  angle. 

The  general  height  level  of  the  vortices  appears  to  be  affected  by 
"bluntnesa"  of  the  nose  and  the  cross  section  shape.  The  vortices  on  the 
V/STOL  nose  (the  "bluntest")  tend  to  lie  vary  close  to  the  surface  and, 
in  fact,  lie  within  the  lateral  confines  of  the  nose  (as  seen  in  a top 
view)  at  10  degrees  of  sideslip.  The  vortices  on  the  F-5  nose  also  tend 
to  lie  close  to  the  top  surface  (within  a body  radius)  to  the  highest 
angle  of  attack  tested,  although  the  leeward  vortex  at  angle  of  sideslip 
is  well  outboard  of  the  flank  of  the  body.  At  the  other  extreme,  the 
elliptical  ogive  with  major  axis  vertical  has  a large  vertical  displace- 
ment of  one  vortex  from  the  top  body  surface,  which  would  prevent  any 
interaction  of  that  vertex  with  a canopy,  for  instance. 

Of  the  three  "research"  shapes  tested,  the  elliptical  ogive  with 
major  axis  horizontal  showed  the  highest  angle  of  attack  for  the  onset 
of  lateral  effects  and  the  least  sensitivity  to  Reynolds  number  effects. 
This  shape  tends  to  fix  the  location  of  separation  more  than  the  others 
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because  of  the  high  qircumferentlal  pressure  gradients  on  the  flanks,  so 
the  variations  in  separation  location  which  cause  lateral  forces  to  occur 
should  be  least  of  the  shapes  tested.  The  forward  portion  of  the  F-5 
nose  has  this  shape,  which  might  be  another  factor  in  its  unusual  force 
characteristios . 

5.2  Vortex  Positions  and  Strengths 

The  tests  of  the  V/STOL  configuration  in  both  the  LRC  V/STOL  tunnel 
and  the  water  tunnel  permit  comparisons  of  vortex  positions  with  the 
theory  described  in  Section  4 and  comparisons  of  the  LRC  results  with 
theory  for  vortex  strength.  The  theory  was  applied  to  the  V/STOL  shape 
using  the  noncircular  body  model  for  35  degrees  angle  of  attack  and  both 
0 and  10  degrees  of  sideslip.  For  zero  sideslip,  an  asymmetric  distur- 
banco  to  the  separation  line  location  was  introduced,  as  discussed  in 
Section  4.2.  The  vortex  center  locations  from  the  LRC  flow  field  data 
were  estimated  by  plotting  radial  and  circumferential  variations  of* 
velocity  components  and  interpolating  for  changes  in  sign  of  tho  cross- 
flow plane  velocity  oompononts. 

The  results  of  the  comparison  for  vortex  positions  are  shown  in 
Figure  19  for  Body  Station  51.9,  which  is  at  the  intersection  of  tho 
strake  leading  edge  with  the  body.  For  the  case  of  zero  sideslip,  the 
agreement  is  reasonably  good.  Both  sets  of  data  and  the  theory  indicate 
little  asymmetry  in  position.  At  10  degrees  sideslip,  the  differences 
are  greater,  with  the  measured  vortex  positions  closer  to  the  body  than 
the  predicted  values.  One  of  the  uncertainties  in  making  this  comparison 
is  the  problem  of  characterd.zing  a distribution  of  vorticity  by  its 
"center" , rather  than  using  the  vorticity  distribution  itself  or  induced 
velocities . 

The  influence  of  Reynolds  number  on  vortex  positions  could  not  be 
determined  with  these  measurements.  Turbulent  separation  was  used  with 
the  theory  of  Section  4,  and  the  predicted  separation  locations  on  the 
body  are  insensitive  to  Reynolds  number.  The  NEAR  wator-tunncl  results 
showed  considerable  variation  in  forces  and  moments  with  a Reynolds 
number  change  from  0.95x10*'  to  S.lxio*",  but  no  vortex  position  data  could 
be  obtained  at  the  higher  tunnel  velocity  due  to  limitations  on  the 
instrumentation.  The  LRC  force  data  show  little  influonco  of  Reynolds 
number  from  2.3x10*’  to  7.4x10*’,  but  no  vortex  position  data  were  obtained 
at  the  lower  Reynolds  number. 
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With  the  velocity  components  known  at  grid  points  such  as  was 

obtained  in  the  LRC  tests,  it  is  possible  to  perform  contour  integrations 
to  obtain  circulation  strengths  for  comparison  with  the  predicted  values. 
A first  approximation  to  such  a calculation  was  made  for  the  case  dis- 
cussed above,  with  10  degrees  of  sideslip,  by  using  two  contours  to 
encompass  the  grid  of  velocity  measurements,  with  the  boundary  between 
the  left  and  right  contours  a straight  radial  line  along  the  plane  con- 
taining the  crossflow  vector  (the  fifth  radial  line  of  points  left  of 
the  vertical  axis  in  Figure  A-17} . From  examination  of  Figure  A-17,  it 
is  clear  that  while  the  outer  part  of  this  common  boundary  is  aligned 
with  the  crossflow  velocity  vectors  and  is  thus  a reasonable  boundary 
between  positive  and  negative  circulation,  the  inner  part  is  not.  As  a 
result,  one  would  expect  the  magnitudes  of  both  the  vortex  strengths  to 
be  somewhat  low  with  this  computation.  The  results  for  vortex  strength 
(r/V)  and  the  predicted  \'aluea  are: 

LRC  Test  Predicted 

right  vortex  +0.520  +0.636 

left  vortex  -0.342  -0.644 

The  predicted  values  show  little  difference  in  magnitude,  whereas  the 
measured  values  show  a considerable  difference,  right  to  left.  A more 
detailed  integration,  particularly  near  the  body,  will  probably  show  an 
increase  in  the  magnitudes  of  the  measured  vortex  strengths,  and  this  i:^ 
planned  in  subsequent  analysis  of  this  data. 

It  is  planned  that  rhe  contour  integration  process  will  be  automated 
so  better  contours  can  be  taken,  and  contours  of  different  size  can  be 
taken  to  determine  the  distribution  of  vortioity. 

5.3  Forces  and  Moments 

It  is  possible  to  compare  predicted  forces  and  moments  on  the  V/STOL 
nose  with  those  measured  in  the  NEAR  water  tunnel  tests.  This  was  done 
for  the  35  degree  angle  of  attack  case  at  both  zero  and  10  degrees  of 
sideslip,  with  the  following  results.  The  predicted  values  were  obtained 
using  a turbulent  separation  criteria  for  separation  location. 
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In 

general,  the  predicted  values  agree  reasonably  well  with  the 

measured 

values  for  the  comparable 

Reynolds  number, 

except  for  pitching 

moment.  The  reason  for  the  disagreement  on  is  not  known,  it  does  not 
have  much  separation-induced  content,  so  the  problem  may  lie  with  the 
force  balance  readings  or  calibration. 


6.  CONCLUDING  REMARKS 

As  part  of  a continuing  study  on  spin  entry  of  high  speed  aircraft/ 
a series  of  analytical  and  experimental  investigations  have  been  performed, 
directed  principally  at  vortex  separation  from  aircraft  nose  shapes  at 
high  angles  of  attack  and  small  sideslip.  The  asymmetric  separation  char- 
acteristics of  the  nose  have  been  identified  as  an  important  factor  in 
the  lateral  aerodynamic  forces  and  moments  that  are  developed  on  many 
fighter-bomber  aircraft  at  incipient  spin  entry  conditions.  Understanding 
and  prediction  of  the  nose  flow  characteristics  will  lead  to  an  ability 
to  calculate  the  lateral  loads  on  complete  configurations. 

A significant  aspect  of  the  investigation  is  the  strong  correlation 
between  the  analytical  work  and  the  execution  of  the  tests  to  obtain  data 
to  check  the  methods  and  provide  understanding  of  the  flow  over  real  nose 
shapes.  In  particular,  it  is  essential  to  obtain  flow  visualization  and 
flow  field  data  under  identical  conditions  as  force  and  moment  data  to 
understand  the  phenomena. 

The  comparisons  between  the  asymmetric  vortex  uhedding  theory  and 
the  various  sets  of  data  have  shown  reasonably  good  correlation  considering 
the  complexity  of  the  flow  phenomena  and  the  fact  that  few  results  are 
available  from  previous  work  on  vortex  shedding  from  noncircular  bodies. 

It  is  necessary  to  do  additional  work  on  the  nose  vortex  shedding 
method  before  It  can  be  incorporated  into  an  overall  computation  method 
for  a complete  aircraft.  The  separation  criteria  need  to  be  examined  in 
more  detail  and  compared  with  additional  data,  particularly  for  three 
dimensional  separation  cases.  The  moans  of  introducing  asymmetric  per- 
turbations for  esses  of  zero  sideslip  at  high  angles  of  attack  need 
additional  examination.  Finally,  there  are  considerably  more  comparisons 
and  evaluations  that  should  bo  done  with  the  LRC  and  water-tunnel  data 
to  evaluate  the  vortex  shedding  analysis. 

When  reasonable  confidence  has  bean  established  in  the  prediction 
method,  it  would  then  be  of  great  interest  to  evaluate  systematically 
aircraft  nose  shapes  to  determine  those  characteristics  which  provide 
the  least  undesirable  or  perhaps  favorable  aerodynamic  behavior.  It 
should  be  possible  to  rationally  design  aerodynamic  characteristics  into 
a fighter-bomber  aircraft  to  provide  desirable  high  angle  handling 
qualities. 
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AFPEKDIX  A 

RESULTS  OF  LANGLEY  RESEARCH  CENTER 
TESTS  IN  VSTOL  TUNNEL 
ON  VSTOL  FIGHTER  MODEL 
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Figure  A. 4.  Yawing  Moment  Variation  with  Angle  of  Attack 
for  Pour  Sideslip  Angles.  Tail  off,  Strakesoff. 
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Figure  A. 12.  Varjation  of  Side  Force,  Yawing  Moment  and  Rolling  Moment  with 
Probe  Position  at  a - 35*,  0=0".  Probe  at  BS  98  and 
Innermost  Radial  Position. 
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Figure  A. 13.  Compariaon  of  Repeat  Probe  Velocity  Data  for  Wing-Body 
at  BS  69,  Bg  ■>  102*,  Rake  Radial  Position  2.54  om  (1  inch) 
out  from  Innermost  Station. 
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Figure  A. 14.  Comperiaon  of  Repaat  Probe  Velocity  Data  for  Wing-Body 
at  BS  98,  Og  ■ 65.5*,  Rake  Radial  Position  33  ran  (13  In.) 
from  Innermost  Position. 
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TABLE  B.  1 


I 

VORTEX  POSITION  DATA  j 

i 

left  vortex  right  vortex  \ 


^v/%q 

*v/^eq 

^v/^eq 

^v/^ea 

I M 

-.083 

.14 

.094 

.14 

1 1 

-.14 

.22 

. 19 

.22 

1 9 

-.20 

.32 

.24 

. 32 

i 1 

-.21 

.40 

.27 

.44 

-.23 

.52 

.27 

.55 

' n 

-.21 

.62 

.25 

.69 

-.  19 

.75 

.24 

.83 

-.21 

.87 

.26 

. 97 

I 9 

-.21 

.94 

.25 

1.06 

j 

-.  23 

1.01 

.23 

1.13 

1 'll 

'■  ™ 

-.089 

.11 

.079 

.11 

-.15 

.23 

.14 

.23 

: 

-.  20 

.35 

.20 

.37 

1 1 

-.24 

.46 

.20 

.57 

-.23 

. 58 

.24 

.63 

' 'i 

-.22 

.68 

. 20 

. 75 

-.20 

.80 

. 18 

. 88 

-.20 

. 89 

. 16 

1.  0 

J 

-.28 

1,00 

. 18 

1.  11 

-.34 

1.06 

.12 

1.27 

-.  15 

.09 

.050 

.09 

1 

-.25 

.19 

. 070 

. 19 

-.32 

.30 

.080 

. 30 

■ -j 

-.38 

.42 

.080 

. 39 

-.44 

.52 

.070 

.48 

,v 

-.46 

.63 

. 060 

. 60 

-.  50 

.74 

.050 

.69 

-.54 

.84 

.020 

. 81 

-.  57 

.86 

n 

. 95 

-.  58 

. 88 

0 

. 98 

! S 

TABLE  B. 1 


Continued 


TABLE  B.  1 

.6 

.7 

.8 

.9 

1.0 

VSTOL  30°  0*^  . 1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 

VSTOL  35°  0°  .1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 

VSTOL*  20°  10°  . 1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 


.57 

.57 

.59 

.59 

. 65 

.65 

' 

.70 

.70 

.74 

r 

.74 

0 

. 18 

no 

.18 

ta 

.38 

da 

ta 

.38 

. 54 

.54 

.61 

.61 

.63 

.63 

.63 

.63 

.64 

.64 

.65 

.65 

f 

.70 

' 

.70 

1 

.76 

1 

.76 

.20 

- 

.20 

23 

.42 

• 

21 

.42 

23 

.56 

21 

.56 

21 

.65 

• 

IS 

.65 

17 

.68 

• 

19 

.68 

14 

. 70 

• 

16 

.70 

16 

.71 

t 

16 

.71 

17 

.75 

t 

12 

.75 

19 

.80 

» 

14 

.80 

.86 

- 

.86 

37 

.24 

• 

056 

.44 

43 

.34 

037 

.56 

48 

.40 

10 

.60 

47 

.45 

11 

.60 

49 

.47 

093 

.61 

49 

.47 

-.084A^0 

.61/.  61 

45 

.45 

-.17/+. 092 

. 65/. 65 

45 

. 5 

-.22/+.  092 

.65/.  65 

45 

. 51 

-/+.092 

.65/.  65 
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TABLE  B.  1 Continued 


VSTOL’ 


VSTOL 


ELLIPTIC 

OGIVE 

MAJOR  AXIS 
HORIZONTAL 


.1 

-.19 

- 

.038 

- 

.2 

-.29 

.43 

.057 

.43 

.3 

-.32 

.54 

.038 

.54 

.4 

-.34 

.63 

-.019/.  15 

.63/.  63 

.5 

-.34 

.65 

-.057/.  17 

.65/.  65 

.6 

-.34 

.66 

-.095/. 15 

. 66/.  66 

.7 

-.34 

.62 

-.095/. 13 

. 69/.  69 

.8 

-.34 

.65 

-/.  13 

. 76/.  76 

.9 

-.34 

.63 

-/.  038 

. 82/. 82 

1.0 

-.34 

.61 

-/o 

. 86/. 86 

.1 

- 

- 

- 

- 

.2 

-.23 

.36 

- 

- 

.3 

-.31 

.48 

.078 

- 

.4 

-.35 

.63 

.078 

.63 

.5 

-.35 

.63 

.088 

.68 

.6 

-.30 

.63 

.088 

.70 

.7 

-.28 

.64 

.078 

.76 

.8 

-.19 

.65 

.097 

.85 

.9 

-.  16 

.62 

.097 

.97 

1.0 

- 

.63 

.097 

1.  12 

.1 

_ 

.12 

- 

.12 

.2 

-.098 

.21 

.098 

.21 

.3 

-.  16 

.26 

.16 

.30 

.4 

-.20 

.35 

.21 

.41 

.5 

-.26 

.41 

.26 

.48 

.6 

-.29 

.47 

.29 

.59 

.7 

-.31 

.52 

.32 

.70 

.8 

-.33 

.59 

.34 

.77 

.9 

-.  31 

.65 

.37 

.83 

1.0 

-.2'J 

.70 

.40 

.87 

*Some  indication  exists  that  a third  vortex  is  forming  under  the  first 
vortex  on  the  right  aide. 
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TABLE  B.  1 Continued 


i 

i 

1 


ELLIPTIC 

o 

o 

0° 

.1 

-.061 

.11 

.061 

.11 

OGIVE 

-.11 

.22 

MAJOR  AXIS 

.2 

.11 

.27 

HORIZONTAL 

.3 

-.15 

.28 

.15 

.39 

.4 

-.20 

.36 

. 20 

.51 

.5 

-.24 

.42 

.26 

.64 

.6 

-.24 

.48 

.29 

.75 

.7 

-.24 

.55 

.34 

.94 

.8 

-.23 

.59 

.40 

1.12 

.9 

-.18 

.64 

.45 

1.  32 

1.0 

-.16 

.70 

.51 

1.57 

ELLIPTIC 

o 

o 

o 

o 

.1 



.15 

.15 

OGIVE 

-.20 

.24 

MAJOR  AXIS 

. 2 

.039 

.24 

HORIZONTAL 

.3 

-.26 

.32 

.039 

.32 

.4  -.33 

.5  -.41 

.6  -.47 

.7  -.57 

.8  -.65 

.9  -.71 

1.0  -.79 


ELLIPTIC  40° 
OGIVE 

MAJOR  AXIS 
HORIZONTAL 


.1 

.2 

.3 


.16 

,22 


TABLE  B, 1 Continued 


ELLIPTIC  30°  0° 

OGIVE 

MAJOR  AXIS 
VERTICAL 


ELLIPTIC  40°  0° 

OGIVE 

MAJOR  AXIS 
VERTICAL 


ELLIPTIC  40°  10° 

OGIVE 

MAJOR  AXIS 
VERTICAL 


.1 

- 

- 

- 

- 

.2 

-.057 

.28 

.048 

.28 

.3 

-.076 

.43 

.076 

. 38 

.4 

-.094 

.60 

.057 

.47 

.5 

-.11 

.7? 

.057 

.53 

.6 

-.13 

.94 

.057 

.59 

.7 

-.15 

1.13 

.094 

.64 

.8 

-.15 

1.34 

.13 

.68 

.9 

-.16 

1.60 

.13 

.74 

1.0 

-.20 

- 

.13 

.77 

.1 

-.060 

.21 

- 

.17 

.2 

-.099 

.44 

.038 

.29 

.3 

-.11 

.74 

.038 

.38 

.4 

-.14 

1.03 

.038 

.50 

.5 

-.18 

1.35 

.049 

.59 

. 6 

22 

1.71 

.058 

.66 

.7 

-.23 

2.02 

.049 

.76 

.8 

- 

2.47 

.049 

.84 

.9 

- 

2.77 

.058 

.90 

1.0 

- 

- 

.077 

.96 

.1 

no 

- 

- 

- 

.2 

data 

.25 

0 

.31 

.3 

. 35 

0 

.49 

.4 

.46 

.04 

.74 

.5 

. 54 

.04 

.99 

.6 

.61 

0 

1.28 

.7 

.63 

-.06 

1.63 

.8 

.71 

-.09 

2.01 

.9 

« 

p 

.79 

-.14 

2.38 

1.0 

. 82 

-.20 

2.75 

% 

ill 

1 


i 


1 


110 


35 


n,  degrees 


Figure  B.l.  Normal  and  Side  Force  Variation  with  Angle  of 
Attack  for  \'/STOL  Nose.  6 = 0. 


Figure  B.2. 


Figure  B.8.  Pitching  and  Yawing  Moment  Variation  with  Sideslip  for 
Elliptical  Ogive  with  Major  Axis  Horizontal. 
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Figure  B.ll.  Normal  and  Side  Force  Variation  with  Sideslip  for 
Elliptical  Ogive  with  Major  Axis  Vertical. 
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Figure  B.14.  Pitching  and  Yawing  Moment  Variation  with  Angle 
of  Attack  for  P~5  Nose.  8-0 
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Figure  B.19.  Force  and  Moment  Changes  on  Elliptic  Ogive  Nose 
Major  Axis  Vertical,  Due  to  Probe  with  and  without  Air. 
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APPENDIX  C I 

MAPPING  OF  CORRBSPOblDING  POINTS  BETVnSEN  | 

j THE  REAL  PLANE  AND  THE  CIRCLE  PLANE  i 

S i 


j 


Introduction 


Calculation  of  the  potential  flow  field  around  an  arbitrary  shaped 
body  and  prediction  of  the  motion  of  free  vortices  in  the  vicinity  of  the 
body  require  that  corresponding  points  in  each  plane  be  readily  obtainable. 
From  Figure  11,  the  position  vectors  in  each  plane  are  written  as 


T + iX 

-ie 


(C.l) 


and 


y + iz 


-ise 


iB 


(C.2) 


Given  the  transformation  from  equation  (1), 


da 

dv  “ 


N - iB_ 

I 

n-0  v" 


(C.3) 


where 


.n 


n ,-in8 
o 


{C.4) 


and  assuming  that  corresponding  points  and  are  known,  the  various 


required  mapping  procedures  are  carried  out  as  follows. 


Map  Known  Points  on  the  Circle  to  the 
Corresponding  Points  on  Body 


Given  the  point  on  the  circle  (Fig.  11),  determine  the  corre- 

sponding point  Cj  on  the  body  as  follows.  From  equations  (C.l)  and 
(C.2) 


dv  - -r^ie"^®de 


(C.5) 


do  ■ -i  da  e^^  + s dp  e^^ 

> (s  dp  COB  p 4-  ds  sin  P) 

+ i(a  dp  sin  p - ds  cos  P) 
■ dy  + i dz 


(C.6) 

(C.6a) 
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I 


I ; 


i 

I ' 


Rewrite  (C.3}  using  (C.5)  and  (C.6} 


! 


da  _ 3 dB  - i da  e^^ 

-i  r e"^®de 
o 


da 

dv 


coB(e  + p)  - ^11  ainte  + p) 


+ i 


sin (6)  + p)  + ^ j|  co8(e  + p) 
o o 


From  (C.3),  define 

R ■ Real  exp 


(f 

' n-0  ^ ‘ 


N 


Therefore, 


I ■>  Imag  exp 


da 


Vn-0 


R - COB  {9  + P)  - ^ f|  8in(9  + p) 

O ■‘’O 


(C.7) 


(C.8) 


(C.9) 


(C.IO) 


(C.ll) 


1 - Bin(0  + p)  + ^ If  CO8(0  + p)  (C.12) 

o o 

Since  R and  I are  functions  only  of  the  Itnown  point  on  the  circle, 
solve  (C.ll)  and  (C.12)  co  get 


da 

r^de 

o 


- R C08(e  4-  P)  + 1 8in(e  + p) 


^ If  - I COB(0  + p)  - R Bin(0  + P) 


(C.13) 


(C.14) 


From  (C.6) , 


t 

I 


I,;'  ' 

% '' 


if)'"' ® + 

■’*  ■ [(^  if)  •^"  ^ - (^) 


(C.15) 


(C.16) 


The  coordinates  on  the  body  are 


Vx  • Vq* 


<C.17) 


•*•  (<J»)ro<Je 


Map  a Specified  Point  Near  the  Body  to  Its 
Corresponding  Point  Near  the  Circle 

Considering  the  following  sketch,  given  a^,  find  Vj 
z 


(C.18) 


s 

a ■ y + iz 


V ■ T + iX 


The  position  of  Oj  uniquely  determines  the  point,  Oj,  on  the  body. 

Prom  the  previous  mapping,  the  corresponding  point  on  the  circle,  Vj, 
can  be  found  either  from  a table  lookup  or  a stepwise  application  of  the 
previous  mapping.  The  procedure  to  locate  Vj  is  to  integrate  along  the 
straight  line  o7a7  (S  • constant)  and  determine  the  corresponding 


curved  line,  Prom  (C.2), 

0 - -is  e^ 


(C.19) 


and  holding  B constant. 


da  ■ -i  ds 
“>  dy  + i dz 


From  (C.l), 


V ■ re 


-10 


(C.20) 

(C.21) 


and 


dv  - -Ir  de  + dr 


- dT  + 1 dX 

Expanding  (C.22)  results  In 

d\J  •*  (dr  cos  6 - r d6  sin  6)  - l(r  d0  oos  6 -f  dr  sin  6) 
Combining  (C.3)  with  (C.9)  and  (C.IO), 

1B„ 


da 

Jv  • 


N 

I ^ 


|^n«0 

From  (C.20)  and  (0.22), 

dv  _ (dr  - Ir  d9)e 
“ -1  ds  0 

" 3i-  H 


« R + II 


-19 


■1(9+0) 


Jr  11  008(9  + 0)  + §1  sln(0  + 0>J 


+ 1 


||  008(9  + 0)  - r ^ Bln(0  + 0) 


(C.22) 


(C.23) 


(C.24) 


(C.25) 


Let 


dv 

3a 


1 

15 

av 


P + IQ  - 


R + II 


(C.26) 
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where 


P - Bin(0  + p)  + r II  co8(0  + p)  {C.27) 

Q - ||  coa{0  + p)  - r II  ■in{e  + p)  (C.28) 

Note  that  F and  Q are  funotlone  of  V/  but  in  each  interval  ds>  v 
can  be  treated  as  constant.  In  the  numerical  integration/  some  error 
is  introduced,  but  the  results  will  converge  as  ds  is  smaller. 

Solving  (0.27)  and  (C.28),  we  get 


dr  _ 
ds 


Combining  (C.22), 


P Bin(d  + P)  Q cos(0  p) 

P oos(0  + p)  - Q 8in(0  + p) 

(C.23),  (C.29),  and  (0.30),  we  get 


(C.29) 

(C.30) 


dT  8in(0  + p)  + Q cos(0  + P)j  cos  0 

“^P  cos(0  + P)  - Q sin(0  + P)j  sin  0^ds' 
dX  sin(0  + p)  - Q cob(0  + p)jsin  0 

-^P  CO8(0  + P)  - Q sin(0  + p)j  cos  0}da 
The  coordinates  of  are 

s 

‘ li 

B 

1 

K " \ f I? 

a„ 


(C.31) 


(C.32) 


(C.33) 


(C.34) 
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1 


From  (C.2} , 


and  from  (C.6}( 


a ■ -is  e 


da  - (a  dS  - i ds)e 


The  transformation,  from  (C.3) 


A - is 

^ .,n 


R -f  il 


where 


v"  . r”a-i"* 


can  aleo  be  written  in  terms  of  (C.36}  and  (C.36) 


(a  dS  - i dsie^^ 

-i.e 


- 8 II  coa(9  + p)  + Il  8in(e  + 0) 

+ i 1^8  II  8in(0  + P)  - If  cos(e  + P) 


From  (C.39)  and  (C.41), 


R ■ 8 H cos (9  + p)  + If  sin(0  + P) 


I - 8 II  Bln(0  + P)  - If  CO8(0  + p) 

where  R and  I are  funations  of  r only  from  (C.39) 
and  (C.43}  results  in 


(C.37) 


(C.38) 


(C.40) 


(C.41) 


(C.42) 


<C.43) 


Solving  (C.42) 


(C.44) 


R C08(e  + $)  + I sin(0  + S) 

R sin(e  + e)  - I co8(6  + 3)  (C.45) 

From  (C.38) 

do  ■ (6  d3  COB  3 -f  d8  sin  3) 

+ l(s  d3  sin  3 - ds  cos  3) 

■ dy  + i dz  (C. 46) 

Therefore, 

dy  >■  8 d3  008  8 + ds  sin  3 (C.47) 

dz  8 d3  sin  3 **  ds  oob  3 (C.48) 

Combining  (C.44)i  (C.45),  (o.47),  and  (C.48) 

dy  "^fR  C0B(e  + 3)  + I 8in(0  + 3^  cos  3 

+ Qr  8in(6  + 3)  “ I 008(0  + 3^  sin  S^dr  (C.49) 

dz  co8(0  + 3)  + I Bin(e  + 3)^  sin  3 

- Qr  8in(e  + 3)  - I 008(0  + 3)3cos  3}  dr  (C.50) 

Finally,  the  coordinates  of  the  point  are 

y*  - yi  j 

which  are  obtained  by  numerical  integration.  Note  that  dy  and  dz 
require  values  of  3 along  the  path  of  integration.  It  has  been  found 
that  in  integration  in  a stepwise  fashion  with  nmell  dr  steps ^ the 


(C.51) 


(C.52) 


® a? 


ds 

a? 
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p from  the  previous  step  can  oe  uaea.  xna-s  introauaes  onxy 
ror  in  the  final  results. 


APPBNOIX  D 

THE  UNSTEADV  PRESSURE  UISTRIBUTIOH  OH 
GROWING  BODIES  UNDER  THE  INFLUENCE  OF 
FREE  VORTICES 


The  Bernoulli  equation,  written  as 

p«  + + p H " p * p H 


(D.l) 


requires  evaluation  of  the  unsteady  pressure  term  at  infinity  and  on  the 
body  surface.  The  velocity  potential  from  equation  (6)  is 


B Real  W(c] 

where  the  complex  potential  from  equation  (17)  is 


(D.2) 


W(o)  - -(V  sin  3) 


i(V  sin  a)  V - 


^ 27  v~-"\r  'liT 


(D.3) 


Now  we  need 


^ ^ [Real  W(o)] 


V ■ T + iX 


(D.4) 


and  consider  one  term  of  (D.3)  at  a time.  Starting  with  the  g-oross- 
flow  term, 


- r®  r*  “ 

W(a)  - -V  sin  6 t[i  + , + ufl  - T-^] 

I'  T+X^  T+X-* 


applying  equaU-on  (D.4)  produces 

r 


■3^  ■ '“'“I 


(D.6) 


(D.7) 


The  a-orossflow  term. in  (D.3)  is 


r 1 "I 

(a)„  ■ -iV  sin  ot  t 1 - — j r + iX  1 + — s j- 

“ T+X'  ' T+X'^ 


(D.8) 


I- 
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which  results  in 


^ . V .in  v4 


The  r term  in  {D.3)  is 
n 


N r 

W(o)j.  - i 

n-1 


- In 


fx  - ^ 

2 J-  ■>  ® J 

n n 

[x  . ,J  + i(i  - A„] 


X - 


n n 


Rewriting  (D.IO)  as 


W(a)„  “ i I ■>5^  [ln{pe^^)  - In  (re  ) 
^ n"l  “ *- 


>] 


where 


T - 


r*  T 
o n 


^ a 


■"n  + ^nj 


X - 


r*  X 
o n 


* X \ * 

n n 


i/a 


q = tan 


- 1 


1 

J 

' r*  X ] 
o n 

n H'' 

r a X 

r x_ 

o n 

^ mm 

2 . 2 

r^ 

+■ 

>- 

mm 

n n 

(x-x„rM-v]‘' 


1/2 


B ■■  tan 


- 1 


X - X. 


T - T. 


Therefore, 


N r. 


Real  W(a)j,  - I 
* n-1 


“ q) 


(0.9) 


(D.IO) 


(D.ll) 


(D.12) 


{D.13) 


(D.14) 


(D.15) 


(D.16) 
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Aesuitiing  is  constant, 


d(t> 


r 

dt 


N r 

y — 

2ir 


n=l 


ds  dc! 
dt  dt 


(D.17) 


which  becomes 


difrr 


dt 


-{t  - T )X  + (X  - X )t 
n n n n 


(X  - XJ* 


{t('/  + X*)  - r*T„}{2X (t_t_  + Xi)  - rfX  - 2r  X_r_} 


on 


n n 


n n 


o n 


o n o 


r / * , 2 » 2 , , 2 

{T(Tn  + X^)  - r^X^}  + 


{X(t'  + x!)  - r'x  }‘ 

n n on 


{X{'t*  + X*)  - rlx  ]{2t{t  -r  + XX)  - r*T_  - 2r  T r } 
n n on n n n n on  one 


n n on  n n on 


(D.iti) 


The  last  term  in  (D.3)  represents  the  expanding  body  radius 


dr 

W(a)^  * (V  cos  a)r^  ln(T  + iX) 


(D.19) 


Writing 


1 + iX  = pe 


iY 


(D.20} 


where 


ln(T  + iX)  = In  p + iY 


(D.21) 


we  get 
dd 


r 

o 

'd't" 


(V  cos  a) In  p 


• ^^o 

^o  ~Sx 


+ r 


o dt 


dr 

o 

dx 


(D.22) 


It  can  be  shown  that  each  component  of  (D.4)  given  by  (D.7), 

(D.9},  (D.IS),  and  (D.2?}  is  finite  at  the  body  boundary,  and  all  but 
(0.22)  disappear  as  v Including  (0.22)  will  not  allow  an  absolute 

pressure  coefficient  on  the  body  to  be  computed.  The  computed  pressure 
distribution  at  any  crossflow  plane  will  be  indeterminate  by  an  dditive 
Constant  factor.  As  done  in  reference  6,  this  factor  can  be  computed 
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by  comparing  the  predicted  pressure  distribution  with  a measured  distri- 
bution and  choosing  a factor  which  produces  good  agreement  between  the 
curves . 

An  alternate  approach  to  this  difficulty  is  to  choose  a different 
singularity  distribution  to  represent  the  growth  of  the  body  radius.  A 
three-dimensional  distribution  of  sources  and  sinks  along  the  body  axis 
to  ^present  the  body  thickness  distribution  is  described  in  reference 
10.  The  procedure  is  to  distribute  souri.’es  and  sinks  along  the  body  axis 
in  such  a manner  that  the  body  .surfa'~e  is  a stream  surface.  This  is 
carried  out  by  matching  body  radius  and  slope  at  a large  number  of  pre- 
scribe'^ points. 

Given  such  a distribution  of  sources,  the  velocity  potential  is 


where 


K 

(X)  = I 


k“l  47t/(x  - Xj^)  * + r^ 


(D.23) 


cos  a Q* 
n 


(D.24) 


reference  length 


= dimensionless  source  strength 


= total  number  of  singularities 


db  (x)  _ V cos  g r ^k  ~ ^ d^ 

dt  ■ 47r  . i r j 2-|  ’/S' 


(x  - Xj^)  + r 


(D.25) 


It  can  be  shown  that 


lim  = 0 


(D.26} 


and  on  the  body  surface, 


db (x)  _ V cos  g y _k 

<4-^  ™ A-rr 


K Qj  (X  - xj^)  + -ad 


' 0"  ‘ 


V 2 n 

X,  ) + r 

k o I 


(D.27) 


Thus,  replacing  (D.22)  with  (D.27)  provides  a method  of  computing  the 
absolute  pressure  coefficient  on  the  body  surface. 

The  complete  unsteady  term  is  the  sum  of  {D.7),  (D.9),  (0.18),  a 
(D.27}.  Since 

^ “ V cos  o„ 


we  can  write 


V cos  a. 


• o dx  ^ o 

^o  “ "dsT  dt  “ ^ “c  "dir 


_ _d_ 

dx  ,, 

dx 

dx 

[dxj 

— ■ = V cos  a, 
dt  c 

Therefore,  on  the  body  surface 


■ -V  sin  6 V cos  -g^j  + V sin  a ^ V cos  -gjS 


N r 

+ 2W 


(t  - T ) 


n'  ■3F 


(X  - 


n'  ~S^ 


(T  - t:„)'  + (X  - 

n n 


■cr  - T r 2Xt  -r-li  + 2Xa 


T r_)  + (Xr 


- 2r  X 

o n dx 


o n dx 


no  n no 


(continued  on  next  page) 


(D.32) 


, , ^ Qfc  (X  - X.  ) + r 

+ v^cos  o I : 


(X  - 


V 2 . * 


Finally,  (0.32)  can  be  simplified  to 

11  - - ^ sin  6 cos  ^ sin  a cos 

VO  o 

, . 2 ? QJR*  ■ 'o  -33^1 

+ 2 cos  a I 3/r 

k-1  ^ V.2I  ' 


4-  2 cos 


i ^ (-■) 


and  can  be  included  in  aquation  (38) . 


(D.33) 


I 


SYMBOIiS 

Fourier  coefficients  of  transformation,  equation  (1) 

normal  force  coefficient  per  unit  length,  equation  (63) 

sideforce  coefficient  per  unit  length,  equation  (67) 

pitching  moment  coefficient,  c ■ pitching  moment/gSj,d  , 
positive  nose  up  ” ° 

yawing  moment  coefficient,  C_  • yawing  moment/qS.  d , 
positive  nose  right 

pressure  coefficient,  equation  (37) 

modified  pressure  coefficient,  equation  (41) 

normal  force  coefficient,  > normal  force/qSj^,  positive 
nose  up  " 

side  force  coefficient,  ■ side  foroe/qS.  , positive  nose 
to  starboard  * “ 

dieuneter 

diameter  of  circle  having  same  area  as  base  of  model 
separation  parameter,  equation  (48) 
diffusion  term,  equation  (58) 
reference  length 

total  number  of  Fourier  coefficients  used  to  describe 
transformation,  also  model  length 

initial  position  of  shed  vortices,  equation  (54) 

initial  image  position  of  shed  vortex,  equation  (55) 

total  number  of  shed  vortices 

local  pressure 

free  stream  static  pressure 

free  stream  dynamic  pressure,  psf 

source  strength 

dimensionless  source  strength 

radial  distance  between  two  points 

radius  of  circle 
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SYMBOLS  (CONTINUED) 


Re^  free  stream  Reynolds  number,  based  on  d^^ 

s radial  dlstauice  in  a-plane,  equation  (A. 2) 

5 reference  area 

Sjg  base  area 

t time 

U local  velocity 

V free  stream  velocity 

v,w  velocity  components  In  real  plane 

v',vf'  velocity  components  in  circle  plane 

W complex  potential 

x,y,z  body  coordinate  system  with  origin  at  the  nose;  x positive 
aft  along  the  model  axis,  y positive  to  starboard,  and 
z positive  up 

X axial  location  of  center  of  pressure 

P 

x„  axial  location  of  center  of  moments 

m 

a angle  of  attack 

Cq  angle  between  free  stream  velocity  vector  and  body  axis 

0 angle  of  sideslip 

V vorticity  flux 

r vortex  strength 

6*  displacement  thickness 

^t  time  increment 

Lx  axial  length  increment 

6 polar  angle  in  v-plane 

V complex  coordinate  in  circle  plane,  t + iX (figure  1) , also 

kinematic  viscosity 

V effective  kinematic  viscosity 

5 boundary  layer  run  length  along  body  surface 

p density 
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SYMBOLS  (CONCLUDED) 


a complex  coordinate  in  real  plane,  y + iz  (figure  li) 

t,X  lateral  and  vertical  coordinates  in  crossflow  plane  of 

transformed  body 

()>  velocity  potential  in  circle  plane,  equation  (5)  , also  bank 

euigle 

♦ velocity  potential  in  real  plane,  equation  (8) 

stream  function  in  circle  plane,  equation  (5) 

4'  stream  function  in  real  plane,  equation  (8) 

Subscripts 

(*)  derivative  with  time 

( } conjugate  of  complex  quantity 

( associated  with  a body  vortex 
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